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SECTION  I 


INTRODUCTION 

Dynamic  si  ability  or  instability  of  clastic  si  rue l urea  has  drawn 
considerable  attention  in  the  past  thirty  years.  The  Importance  of  the 
subject  lies  primarily  in  the  constant  demand  for  lightweight  efficient 
structures,  which  requires  a  good  understanding  of  how  structures  respond 
to  loads  that  induce  dynamic  effects. 

The  term  "Dynamic  Stability"  encompasses  many  classes  of  problems 
and  it  has  been  used,  by  the  various  investigators,  in  connection  with  a 
particular  study.  Therefore  it  is  not  surprising  that  there  are  various 
interpretations  of  the  meaning  of  the  term. 

The  class  of  problems  falling  in  the  category  of  parametric  excitation 
are  the  best  defined,  conceived  and  understood  problems  of  dynamic  stability. 

An  excellent  treatment  and  bibliography  can  he  found  in  the  book  by  V.  V. 
Bolotin  [1].  As  a  matter  of  fact,  Bolotin  applies  the  term  "Dynamic  Stability" 
only  to  problems  of  parametric  excitation. 

In  general,  problems  which  deal  with  stability  of  motion  have  concerned 
for  many  years  researchers  in  many  fields  of  engineering.  Definitions  for 
stability  and  related  criteria,  as  developed  through  the  years,  are  given  by 
T.  J,  Stoker  [2].  Additional  references  [?.-6]  provide  the  necessary  details 
as  applied  to  a  variety  of  problems.  Some  of  these  criteria  find  wide  uses 
in  problems  of  stability  and  control  of  aircraft  [7],  control  theory  [8],  and 
many  problems  in  fluid  mechanics,  combustion  and,  to  a  certain  extent,  elastic 
structures  under  dynamic  conditions. 

Moreover,  many  authors  refer  to  problems  of  the  "follower  force"  type 
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as  problems  of  dynamic  stability  [9-10].  The  primary  reason  for  this 
Is  that  critical  conditions  can  be  obtained  (In  many  cases)  only  through 
the  use  of  the  "kinetic"  or  "dynamic"  approach  to  static  stability  prob¬ 
lems  (flutter  instead  of  divergence  type  of  instability). 

In  addition,  problems  of  aeroclasiic  instability  and  flow-induced 
instability  (fluid  flowing  through  pipes)  also  fall  under  the  general 
heading  of  dynamic  stability. 

A  large  class  of  structural  problems  that  has  received  attention 
recently  and  does  qualify  as  a  category  of  dynamic  stability  is  that  of 
impulsively  loaded  configurations  and  configurations  which  are  suddenly 
loaded  with  loads  of  constant  magnitude  and  infinite  duration.  These 
two  types  of  loads  may  be  thought  of  as  mathematical  idealizations  of 
blast  loads  of  (a)  large  decay  rates  and  small  decay  times  and  (b)  small 
decay  rates  and  large  decay  times  respectively.  In  addition,  these  types 
of  loads  may  be  considered  as  representative  of  impact  loads.  Of  course, 
there  are  other  physical  explanations  for  such  loads.  For  both  types  of 
loads  mentioned  above  the  concept  of  dynamic  stability  is  related  with 
the  observation  that  for  sufficiently  small  values  of  the  loading  the 
system  simply  oscillates  (linearly  or  non-linearly)  about  the  near  static 
equilibrium  point  and  the  corresponding  amplitudes  of  oscillation  are 
sufficiently  small.  If  the  loading  is  increased,  some  systems  will  experi¬ 
ence  large  amplitude  oscillations  or  in  general  divergent  type  of  motion. 

For  this  phenomenon  to  happen  the  configuration  (turns  out)  must  possess 
two  or  more  static  equilibrium  positions  and  "tunnel ing-through"  [11] 
occurs  by  having  trajectories  that  can  pass  through  an  unstable  static  equi¬ 
librium  point.  Consequently,  the  methodologies  developed  by  the  various 
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investigators  (including  criteria  and  estimates)  are  for  structural 
configurations  that  exhibit  snap- through  buckling,  when  loaded  quasi- 
atatically. 

Solutions  to  such  problems  started  appearing  in  the  open  literature 
in  the  early  1950's.  Hoff  and  Bruce  [12]  considered  the  dynamic  stability 
of  a  pinned  half-sine  arch  under  a  half-sine  distributed  load.  The  ideal 
impulse  problem  as  well  as  the  case  of  a  suddenly  applied  load  with  constant 
magnitude  and  infinite  duration  were  considered  in  this  report.  Budiansky 
and  Roth  [13]  in  studying  the  axisymmetric  behavior  of  a  shallow  spherical 
cap  under  suddenly  applied  loads  defined  the  load  to  be  critical,  when  the 
transient  response  increases  suddenly  with  very  little  increase  in  the 
magnitude  of  the  load.  This  concept  was  adopted  by  numerous  Investigators 
(see  [14])  in  the  subsequent  years  because  it  is  tractable  to  computer  solu¬ 
tions.  This  concept  was  used  by  Budiansky  and  Hutchinson  [15]  in  estimating 
the  critical  load  (suddenly  applied)  for  imperfection  sensitive  systems. 

Through  this  criterion,  they  related  the  dynamic  critical  load  to  the  static 
one.  Finally,  the  concept  was  generalized  in  a  subsequent  paper  by  Budiansky 
[16  ]  in  attempting  to  predict  critical  conditions  for  imperfection  sensitive 
structures  under  time  dependent  loads.  Conceptually,  one  of  the  best  efforts 
in  the  area  of  dynamic  buckling,  under  impulsive  and  suddenly  applied  loads, 
is  the  work  of  Hsu  and  his  collaburators  [17-20].  In  his  studies,  he  defined 
sufficiency  conditions  for  stability  and  sufficiency  conditions  for  instability, 
thus  finding  upper  and  lower  bounds  for  the  critical  impulse  or  critical 
sudden  load.  Independently,  Simitses  [21]  in  dealing  with  the  dynamic  buckling 
of  shallow  arches  and  spherical  caps  termed  the  lower  bound  as  a  minimum  possible 
critical  load  (M^CL)  and  the  upper  bound  as  a  minimum  guaranteed  critical  load 
(MGCL). 
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The  purpose  of  the  present  work  is  to  investigate  the  concept  of 
dynamic  stability  of  structural  elements  subjected  to  step-loads  and 
develop  the  related  criteria  and  estimates  for  finding  critical  conditions. 
The  step  load  consists  of  a  suddenly  applied  load  of  constant  magnitude 
and  finite  duration  tQ,  and  the  investigation  will  include  the  two  ex¬ 
treme  cases  of  t  -•  00  and  t  -*  0  (ideal  impulse).  Moreover,  the  effect  of 
00 

various  parameters  (small  damping  and  preloading)  on  the  critical  con¬ 
ditions  is  studied.  The  developed  solution  methodology  is  demonstrated 
through  simple  mechanical  models  of  one-and  two-degrees  of  freedom.  Section 
I  deals  with  the  introduction  to  the  subject  of  dynamic  stability  of  struc¬ 
tural  elements.  Section  II  presents  a  detailed  static  buckling  analysis  of 
the  models.  Section  III  presents  the  concept  of  dynamic  stability  for  im¬ 
pulsively  loaded  systems,  which  is  demonstrated  through  the  same  models. 
Section  IV  gives  a  detailed  treatment  of  systems  loaded  suddenly  with  con¬ 
stant  load  of  infinite  duration.  Section  V  extends  these  concepts  for  the 
case  of  suddenly  applied  loads  of  constant  magnitude  and  finite  duration. 
Sections  VI  and  VII  deal  with  the  study  of  the  effects  of  preloading  and 
small  damping  on  the  dynamic  critical  conditions. 


SECTION  II 


where  is  the  total  potential,  and  the  character  of  equilibrium  (stable 
or  unstable)  by  the  sign  of  the  second  derivative. 


Fig.  2-1  Geometry  and  Sign  Convention  for  Model  A. 
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The  static  equilibrium  points  are  characterized  by 


p  *  [./l+sine  -  ./U-sinQl  for  6  i  0  (2.3) 

Vl+sinG 

Note  that,  for  9^  *  0  equilibrium  is  characterized  by 

either  9=0 

or  p  =  cot©  („/l+sin6  -  l)/^/l+sinQ  (2.4) 

Equilibrium  positions  are  plotted  on  Fig.  2  as  p  versus  0-6o  for  various 
values  of  the  geometric  imperfection  The  stability  test  reveals  that 

the  dashed  line  positions  are  stable,  while  the  solid  line  positions  are 
unstable  and  snapping  (violent  buckling)  takes  place  through  the  existence  of 
a  limit  point.  Also  note  that  positions  characterized  by  negative  values 
for  8t8o  (not  shown  herein)  are  stable  and  there  is  no  possibility  of  buck¬ 
ling.  Therefore,  our  interest  lies  in  the  area  of  8  >0  and  0-0  >  0. 

o  o 

It  is  assumed,  for  simplicity,  that  the  bars  are  weightless  and  the  mass  m 
of  the  system  is  concentrated  on  the  joint  1} . 


Model  B.  Load  Eccentricity  Model 

Model  B,  shown  on  Fig.  3,  is  represent aL ive  of  eccentrically  loaded 
structural  systems,  exhibiting  limit  point  instability.  The  bar  is  rigid 
and  of  length  L;  the  spring  is  linear  of  stiffness  k  and  the  load  eccentricity 
is  denoted  by  e.  The  static  stability  analysis  for  this  model  is  presented  in 
Ref  [23].  The  bar  is  assumed  to  be  weightless  and  the  mass,  m,  of  the  system 
is  concentrated  on  the  top  of  the  rod,  point  I). 
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Fig.  (2.3)  CccMiii'try  mid  Sfgn  Convention  for  Model  B 


The  expression  for  (lie  total  potential  is  given  by 


U, ^  »  ~ka^sin^9  *  1*1,  (I  *  cos 9  +  ~  sin  8) 

I  Z  1 4 


In  nondimens  iona  1  ir.ed  form  tliis  expression  becomes 

"r  2 

II, r  ^  *  sin  9  -  2p  (I  *  cos  0  +  e  sin  8) 

jkn 


(2.5) 


where  p 


PI. 

ka 


is  the  iioml inieiiR  iona!  Izcd  load  and  c  =  —  is  the  nondlmen* 

la 


slonall/.ed  eccentricity. 


Furthermore,  equilibrium  Is  characterized  by 


=  sinO  (2.61 

p  LanO+G 

and  critical  load  (see  Ref.  [23]l  is  given  by 

Pcr  =  [U'e  2/3l'3/2  (2.7) 

2.4  Model  C.  A  Snap-Through  Model 

Consider  the  model  shown  on  Fig.  (2-4>‘,  which  consists  of  three  equal 
length  rigid  bars.  The  three  bars  are  pinned  to  each  other,  and  they  are 
connected  with  rotational  springs  of  stiffness  P  (linear).  The  left  bar 
is  pinned  onto  an  immovable  support,  while  the  right  bar  is  pinned  onto  a 
movable  support,  which  in  turn  is  connected  to  a  wall  through  a  linear 
extensional  spring  (horizontal)  of  stiffness  k  (linear).  The  middle  bar 
is  originally  horizontal  and  the  loading  consists  of  two  equal  concentrated 
forces,  P,  applied  at  the  ends  of  the  middle  bar  and  remaining  vertical. 

The  original  angle  between  the  horizontal  line,  joining  the  supports  and  the 
end  bars  is  a.  The  angle  between  the  horizontal  and  the  left  bar  in  a 


deformed  state  is  0,  while  the  angle  between  the  horizontal  and  the  right 
bar  Is  <p. 
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This  is  a  two-degree-of- freedom  model,  and  as  it  will  he  seen,  from 


the  ensuing  discussion,  certain  new  features  enter  into  the  solution. 

Assuming  that  a,  0,  and  are  small  angles  such  that  their  sine  can 
be  well  approximated  by  the  angle  itself  and  the  cosine  by  one  minus  half 
of  the  angle  souared,  then  the  expression  for  the  total  potential  is 
approximated  bv 


U*  -  P  [5G2  +  Sep 2  -  2a2  -  »p0  -  20a  -  2pa] 

2 

+  Y  (a2  -  o2  -  92  +  9O)2  -  PL  (2a  -0  -  tp) 


(2.6) 


Note  that  for  the  above  equation  to  be  dimcnsonally  correct  the  following 

order  of  magnitude  must  apply  to  the  various  terms.  If  a,  0  and  c p  are 

2  2 

taken  to  be  small  numbers  fof  order  0) ,  then  P  must  he  of  order  6  x(kL  ), 
PL  of  order  6^x  (kl.2l  and  l)^  of  order  6^c(kL2) . 

The  independent  variables  are  y  and  G  and  the  symmetric  response  mode 
is  characterized  by  Q  =  cp  . 

New  variables  r  and  s,  are  introduced  such  that  the  symmetric  response 
is  characterized  by  s  *  0. 

These  are 


0 


V 


(r+s) 


(2,91 


or 


r 


(e  +  V) 


s  -  ((p  -  0) 


(2.10) 


2 

where  P  ■  P/kL  is  a  nondimens ional ized  rotational  spring  stiffness  para¬ 
meter.  Through  introduction  of  additional  nondimenslonalized  parameters 


U_ 


P2k!.2 


(2.11'' 
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and  by  letting 


a  -  dA.  (2.12) 

the  following  expression  for  the  total  potential  is  obtained. 

u£  -  (r2  +  9s2  -  1JK  r  +  A)  +  -  (A  -  r2-  3s2)2  -  2p  (JR  -  r)  (2.13) 

Note  that,  on  the  basis  of  the  nondimens  ionalization,  the  order  of 

2 

magnitude  of  the  new  parameters  is  as  follows:  (a)  0  is  of  order  6  ;  (b) 

UT  is  of  order  one;  (c)  r  and  s  are  of  order  one;  (d)  A  is  of  order  one;  and 
(e)  p  is  of  order  one. 

It  can  easily  be  shown  that,  if  B  =  0  (no  rotational  springs),  the 
system  is  unstable  for  zero  load,  P  *  0,  and  thus  the  0  *  0  case  is  ex¬ 
cluded  from  the  present  discussion,  which  also  allows  the  nondimens ionaliza¬ 
tion  given  by  Eqs .  (2-12)  and  (2-13)  (division  by  a  nonzero  number).  In 
the  case  0  *■  0.  upon  t lie  application  of  a  very  small  load  P,  the  system 
snaps  through  an  asymmetric  mode  (s  ^0). 


By  introducing  a  new  load  parameter 

Q  =*  P  -  Jh 


(2.16) 
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the  equilibrium  equations,  Eqs.  (2 - IA )  and  (2-lr>)  become 


(A  -  l  -  r2  -  3s2)  r  =  Q 

s(A  -  3  -  r2  -  3s2)  =  0  (2.17) 

There  are  two  possible  solutions  to  Eqs.  (2.14) 

(i)  Symmetric  response  s  =  0,  and 

(A  -  1  -  r2)  r  -  Q  (2.18) 

(ii)  Existence  of  asymmetric  response,  s  ^  0 

2  2 

A  ■  3  *  r  +  3s 

2r  -  Q  (2.19) 


The  equilibrium  positions,  Eqs.  (2-17)  are  plotted  on  Fig.  2-5  for 
all  Q,  and  on  Fig.  2.6  as  a  load-def lection ,  Q-r  curve. 

On  the  basis  of  the  results  and  by  performing  the  stability  test 
(second  derivatives),  the  following  conclusions  are  drawn. 

(a)  For  A  <  1  there  is  no  possibility  of  buckling  and  asynmetric 
response  (s  i  0). 

(b)  For  l  £  A  £  3  the  response  is  symmetric  (s  a  0)  and  buckling 
occurs  through  the  limit  point  (pt.  C  on  Fig.  2-6),  positions  between  A  and 
C  are  stable,  and  the  critical  load  is  given  by 


pcr’/A  +  2  ^ 


(2.20) 
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(c)  For  3  <  A  <  4  there  is  a  possibility  of  asymmetric  modes,  but 
for  this  range  of  A-values,  point  B  is  to  the  left  of  point  C  (Fig.  2-6), 
i.e. 

ft  ■'  1  >  V  A  -  3  (2.21) 

and  buckling  still  occurs  through  the  limit  point.  Therefore,  P  is  given 
by  Eq.  (2.20). 
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(d)  For  A  >  4  buckling  occurs  through  the  existence  of  unstable 
bifurcation  (Pi.  It),  positions  A  to  IS  are  stable,  positions  1SC0  and  HO  are 


unstable,  and  the  critical  load  is  given  by 


Qcr  -  lfl\  -  3  ,  or  (2-22) 

P  -  Ji\  +  2  J\  -  3 
cr 


Pig.  2.6  Load-Displacement  Curve  for  Model  C. 
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SECTION  III 


IDEAL  IMPULSE 


Statement  of  the  Problem: 

Given  a  structural  configuration  subjected  to  loads  that  fall  under 
the  category  of  an  ideal  impulse,  is  there  a  possibility  of  dynamic  in¬ 
stability,  and  if  yes,  what  is  the  level  of  the  impulse  for  this  occurrence? 
Clearly,  before  one  can  estimate  the  critical  condition,  one  must  define 
what  is  meant  by  dynamic  instability.  Some  systems,  regardless  of  the 
magnitude  of  the  impulse,  will  simply  oscillate,  linearly  or  nonlinearly, 
about. the  null  positions.  Others,  if  the  impulse  is  small  will  oscillate, 
but  as  the  impulse  increases  they  have  a  tendency  to  move  away  from  the  null 
position.  Finally,  there  is  a  group  of  systems,  which  behave  a  little 
differently  from  the  second  group.  For  small  impulses  they  simply  oscillate, 
but  as  the  impulse  increases  a  certain  value  is  reached  at  which  and  beyond 
which  the  system  will  oscillate  but  with  a  higher  amplitude  than  that  for 
impulses  smaller  than  this  certain  value.  Kegardless,  though,  of  the  par¬ 
ticular  ca *e  one  can  always  associate  the  magnitude  of  oscillation  amplitude 
to  the  value  of  the  impulse.  In  such  cases  a  critical  impulse  can  be  related 
to  a  maximum  amplitude  (this  is  the  case  in  deflection  limited  designs). 

There  is  always  a  relation  between  the  system  behavior  under  impulsive 
lo*ds  and  the  static  equilibrium  characteristics  of  the  system.  In  the 
groups  discussed  above,  the  first  corresponds  to  systems  which  have  unique 
and  stable  load -deflection  equilibrium  states,  or  a  stable  post-buckling 
bifurcational  branch.  The  other  two  correspond  to  systems  that  exhibit 
either  limit  point  instability  or  an  unstable  bifurcational  branch  in  the  post- 
buckling  region.  The  emphasis,  therefore,  is  placed  on  the  last  two  groups 
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and  criteria  for  stability  and  estimates  of  critical  conditions  are 
presented  through  the  models  discussed  in  Section  II. 


Model  A 

Consider  Model  A  shown  on  Fig.  2-1,  assume  that  the  bars  are  weight¬ 
less  and  that  the  mass  of  the  system  is  concentrated  at  point  B.  Further¬ 
more,  assume  that  the  load,  P,  is  suddenly  applied  with  very  short 
duration  time,  T  ,  and  that  the  impulse  (PTq),  is  imparted  instan¬ 
taneously  into  the  system  as  initial  kinetic  energy. 

Through  impulse -momentum  one  obtains  the  following  relation 


<PT0> 


(3.1) 


where  0  *  dC/dt  at  0=8. 

o  o 

Since  the  system  is  conservative,  then 


+  T°  =  cons  t .  =  l” 


(3.2) 


where  denotes  the  total  potential  "under  zero  load",  T°  is  the  kinetic 
energy ,  given  by 

T°=jmL202  (3.3) 

Note  that  T°  is  the  Initial  kinetic  energy  imparted  instantaneously 
by  the  impulsive  load. 

Finally  the  expression  for  is  given  by 

U°  -  [v'l  +  sln0  -  Jl  +  sin0o]2  (3.4) 


16 


In  order  to  understand  the  "concept '*  of  stability  or  instability  under 
an  impulsive  load  for  this  one  decree  of  freedom  model,  let  us  consider 
Fig.  3.1,  which  is  a  plot  of  versus  0-6^. 

According  to  Eq .  (3.2),  and  since  T°  is  positive  definite,  motion 
is  possible  if  and  only  if 


T 


o 

i 


(3.5) 


This  implies  that,  for  a  given  initial  kinetic  energy,  Eq.  (3.3),  and 

consequently  a  given  impulse,  say  =  D  (see  Fig.  3.1  -  total  potential 

presented  in  nondimens ional ized  form),  motion  is  confined  in  the  region 

0  T<  0  _  0  <  0TT-  It  is  clearly  seen  then  that,  as  long  as  T?  =  D  < 

X  p  II  i 

fj2  -  +  ainG^)^,  the  motion  of  the  system  is  bounded  and  it  contains 

only  the  stable  "zero  load"  static  equilibrium  point,  ».  Such  a  motion 
is  termed  "unbuckled".  For  the  motion  to  cease  to  be  "unbuckled",  i.e. 
to  become  unbounded  and  cease  to  include  only  the  inil ial  stable  static 
equilibrium  point,  It,  D  must  be,  at  least,  equal  to  t lie  value  of  0°  at 
the  unstable  point  C.  Then  that  point  (C)  can  be  reached  with  zero 
velocity  and  the  motion  can  become  unbounded.  Clearly  if  D  is  even 
slightly  higher  than  the  0°  -  value  at  point  C,  the  motion  does  become 
unbounded  and  it  can  contain  other  static  equilibrium  points,  such  as 
point  C.  Such  a  motion  is  called  "buckled"  and  a  critical  condition 
exists  when  the  impulse  is  large  enough  to  satisfy  the  relation 

T^-U^(C)  (3.6) 

Introducing  nondimens ional ized  time  and  load  parameters 
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(3.8) 


(3.9) 


Use  of  Eqs.  (3.6),  (3. A)  (3.3)  and  (3.9)  yields  the  expression  for 
the  critical  impulse  or 

(PTo^cr  *  2UA  ‘  +  »in0o]  /  si"  9q  (3.10) 

Two  observations  are  worth  mentioning  at  this  point:  (a)  Because 
this  is  a  one  degree  of  freedom  model,  the  critical  impulse  (PT0)C  given  by 
Eq.  3.1  represents  both  the  minimum  possible  (MPCL)  and  minimum  guaranteed 
(MGCL)  values  as  defined  in  Ref.  21.  In  other  words  Hsu's  (17-18)  "sufficient 
condition  for  stability  and  instability"  boundaries  are  coincident;  (b) 
Although  the  concept  presented  so  far  is  clear  and  it  leads  to  a  criterion 
and  estimate  of  the  critical  condition,  it  might  be  impractical  when  applied 
to  real  structures.  In  the  particular  example?  shown  so  far,  it  is  clear  that, 
according  to  the  presented  concept  of  dynamic  instability,  "buckled"  motion 
is  possible  if  the  system  is  allowed  to  reach  the  position  8  ■  n/2.  In  many 
cases  such  positions  may  be  considered  excessive,  especially  in  deflection- 
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Potential  versus  6-9  (Model  A). 


limited  designs.  In  such  cases,  if  0  cannot  he  larger  than  a  specified 

value,  then  the  allowable  impulse  is  smaller  and  its  value  can  be  found 

from  F.q.  (3.6),  if  C  is  replaced  by  the  maximum  allowable  value  of  6,  say 

9  In  this  case 
u 

(PT  )  -2[/l  +  s  in0  -  v‘l  +  sin  0  ]/sin  9  (3.11) 

ocr  L  o  o 


Related  with  this  discussion  is  the  broad  definition  of  stability 
proposed  by  Hoff  (22):  "A  structure  is  in  stable  state,  if  admissible 
finite  disturbances  of  its  initial  state  of  static  or  dynamic  equilibrium, 
are  followed  by  displacements  whose  magnitudes  remain  within  allowable 
bounds  during  the  required  lifetime  of  the  structure." 


Model  B 

Following  the  same  notion  as  for  Model  A,  the  expression  for  the 
critical  impulse,  (pt  )  given  by 


(PT  ) 
r  o  cr 


L 

e 


(3.12) 


where 


T 


■(*?-) 


2.1/2 


(3.13) 


and  l  is  the  bar  moment  of  inertia  about  the  hinge. 

The  expression  for  t he  total  potential,  the  kinetic  energy,  and 
impulse-momentum  are  given  by  (respectively) 


U„,  =  4  ka^  sin^9  -  PL(1  *  cos  0+7*  sin  0) 

I  L  »< 


(PT  ) 
o 


i,  dQ 
e)  dt 


(3.14) 

(3.15) 

(3.16) 
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where  I  le  the  moment  of  inertia  of  the  mass  ra  about  the  hinge. 


In  this  case  also,  if  the  design  is  deflection  limited,  ©  ,  then 

la 

~  ka  musi  be  replaced  by  U°  (©)  =  j  ka^Bin2®^  and  the  allowable 
impulse  in  this  case  is  given  by 


^To^al  lovable 


“  sin©, 
e  I, 


(3.17) 


Another  similar  approach  that  may  be  used  for  finding  critical  impulses 
is  presented  herein  through  Model  B.  It  makes  use  of  Eq.  (3.2),  but  in¬ 
stead  of  associating  critical  conditions  with  characteristics  on  the 
"zero  load"  total  potential,  the  critical  condition  is  associated  with 
characteristics  of  the  phase  plane. 

Eq.  (3.2)  for  this  model,  in  nondimens ionalized  form  becomes 

02  +  8  in2  9  =  [<P\>^(2  (3.18) 

where  0  »  de/dr. 

Clearly,  if  £(PTQ  )^J2  <  1»  (3.18)  denotes  a  closed  curve  about  the 

null  position  Q  ■  0  *  0)  in  the  phase  plane.  In  this  case  the  motion  is 
called  "unbuckled".  (See  Fig.  3.2).  When  ^ptq  =  1,  Eq.  (3.18)  denotes 

a  curve  which  can  escape  the  closed  loop  and  thus  the  motion  becomes  "buckled". 
Therefore, 


(PT  ) 
o  cr 


(3.19) 


Model  C 

In  order  to  simplify  matters,  it  is  assumed  that  the  three  bars  are 
weightless  and  that  the  Impulse  is  imparted  into  the  system  through  two 
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Fig.  3.  2.  Phase  Plane  Curves  for  Model  B 

masses,  m,  at  points  B  and  C.  This  is  also  used  in  determining  the  initial, 
kinetic  energy  of  the  system.  Thus 

(3-20) 

and 

Ti  "  (Jf)|r  •  S 

where  T  is  the  small  duration  time. 


Introduction  of  the  new  nondimens ionalized  parameters 


T  * 


and  -  T  /  , 


(3.21) 


yields 


(3.22) 


The  concept  of  dynamic  stability  is  similar  to  the  one  used  for  Models  A 
and  B. 

0°  +  T°  =  T°  (3.23) 


A  critical  condition  exists,  if  the  impulse,  (pt^) ,  imparts  sufficient 

kinetic  energy  into  the  system  so  it  can  reach,  with  zero  velocity,  an 

unstable  static  equilibrium  point  on  the  "zero-load,,  total  potential, 

0°  and  thus  the  motion  can  become  unbounded. 

T 

Clearly,  then  from  Eqs.  (3.22)  and  (3.23) 

1/2 

(pT0)cr  ‘  C~r°)  =  (unst.  st.  pt.)  (3.24) 


Thus,  before  (pt  )  can  be  found,  one  must  have  knowledge  of  all 
stationary  points  and  of  their  character  (stable,  unstable,  etc.).  Starting 
with  the  expression  for  U°  [Eq.  (2.11)  with  p  =  0],  requiring  equilibrium, 
and  performing  a  stability  analysis  the  following  results  are  obtained. 

Stationary  nolms  and  their  charnel  or 


Pt.  I  at  (V'A,  0) 

Pt.  2  at  [(-v^A  ♦  VA-4)/2.0] 

Pt.  3  at  [f -A-  ^4)/2,0] 

Pt.  4  at  i-jKTi  ,  4/2~> 

Pt.  5  at  l-Jfjl  .  -Jk  -  4/2] 


Stable  (Relative  Minimum) 
Unstable  (Relative  Maximum) 
Stable  (Relative  Minimum) 
Unstable  (Saddle  Point) 
Unstable  (Saddle  Point) 
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c 


* 


t 


I 


We  observe  that  there  are  two  possibilities  of  the  motion  becoming  un¬ 
bounded;  (a)  by  reaching  Pt.  2  with  aero  velocity,  in  which  case  the  system 
will  definitely  move  toward  Pt .  3  (far  stable  point); in  this  case  the  corre¬ 
sponding  critical  impulse  is  termed  (MGCL)  because  U°  (Pt.  2)  is  the  Largest 
of  all  U°  (Stationary  pts.) .  Thus,  regardless  of  the  path  of  motion  in  the  r, 
a  space,  starting  from  the  motion  will  enclose  at  least  one  unstable 

point  and  will  become  unbounded;  (b)  by  reaching  either  Pt.  4  or  Pt.  5,  in 
which  case  there  is  a  single  possibility  of  the  motion  becoming  unbounded  and 
enclosing  one  unstable  point;  note  that,  in  this  case,  there  is  a  possibility 
of  unbounded  motion,  but  this  is  not  guaranteed,  because  the  system  can  simply 
oscillate  in  the  r,  s,  space  bounded  by  lines  of  equal  potential,  to  U°  (Pt.  4) 
enclosing  only  the  near  static  stable  point,  Pt .  1,  and  either  of  the  "saddle" 
points,  Pt.  4  or  Pt .  5;  thus,  the  corresponding  critical  impulse  is  called 
(MPCL) . 

(MPCL)  and  (MGCL)  denote  upper  and  lower  bounds  of  (pt  ) 

o  cr 

By  Eq.  (3.24), 

(MPCL)  (PT0)cr  *  3<A  *  1  (3.25) 

(MGCL)  (PTo)cr  "  [5A-2-3^2-4A  +  £  (A+2 +Jfr2  -  4A)2]1/2 

Note  that  critical  conditions  can  exist  only  for  A  ^  4.  If  A  i  4,  there 


is  no  far  stable  static  equil.  pt .  on  the  "zero-load"  total  potential. 
The  results  are  plotted  on  Fig.  3.3. 
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SECTION  IV 

CONSTANT  LOAD  OF  INFINITE  DURATION 

.The  concept  of  stability,  for  this  type  of  loads  is  also  based  on  the 
definition  of  "buckled"  and  "unbuckled"  motion  as  for  the  ideal  impulse 
Moreover,  the  criterion  for  stability  and  the  related  procedure  for  estimat¬ 
ing  critical  conditions  will  be  presented  through  the  same  three  models. 

Model  A 

For  this  case,  the  sum  of  the  total  potential  and  kinetic  energy  is  zero 


UTP  +  TP  «  0  (4.1) 

-P 

Fig.  4.  1  shows  plots  of  UT  versus  0-9q  for  various  values  of  the  applied  load, 
p.  It  is  seen  from  this  figure  that  for  p  <  0.432  motion  ia  confined  between 
the  origin  and  9-0o  *  A,  or  the  motion  is  "unbuckled".  A  critical  condition 
exists  when  the  motion  can  become  unbounded  by  including  position  A"  (buckled 
motion).  Thus,  the  critical  load  is  found  by  solving  the  following  equations 

uj?  —  0  and  p  -  LrA  +  sin  9  -  ,/l  +  sin  0  ]  - — ■ —  (4.2) 

yr+  sin  e 


with 


,2-P 

duT 


d0 


<  0 


-P 

The  inequality  condition  ensures  that  *  0  at  an  unstable  equilibrium 
point.  The  simultaneous  solution  of  Eqs.  (4.2)  yields  the  dynamic  critical 
load  and  the  corresponding  position  of  the  unstable  static  point,  A"  on  Fig. 
4.1. 
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The  results,  for  this  model,  are  presented  graphically  on  Fig.  4.2. 

For  this  load  case,  also  (MPCL)  and  (MCCT.)  are  one  and  the  same  level. 


Mpdel.B 

For  this  model,  the  critical  condition  is  estimated  through  the  simul¬ 
taneous  solution  of  the  following  equations 

ull  *  -ka^sin^  9  -  PL(1  -  cos  0  +  —  sin  9=0 
T  Z  L 

ka^  sin  0  cos  9  -  PL  (sin  9  cos  0)  =  0  (4.3) 

Li 

2  P 

and  — ■■  <  0  (at  the  solution) 
d0 

The  results  are  presented  graphically  on  Fig.  4.3  as  p  versus  e/L. 

cr 

For  this  model,  critical  conditions  are  also  related  to  characteristics  of 
the  phase  plane  curves.  Fig.  4.4  shows  plots  of  Ur),  versus  0  for  various 
values  of  p,  as  well  as  6  versus  6  (phase  plane). 

The  energy  balance  equation,  in  nondimens ional i zed  form  is 

•2  2  e 

9  +  sin  9  -  p  (1  -  cos0  +  —  sin0)  =  0  (4.4) 

li 

Clearly  from  Fig.  4.4.  if  p  <  p  the  motion  is  "unbuckled",  which 

cr  oo 

implies  that  for  these  loads  the  system  simply  oscillates  about  the  static 
stable  equilibrium  position  (A^  or  •  The  load  p  is  critical,  (p  )  if 

croo 

the  motion  can  escape  through  point  ('buckled  motion") . 

Note  that  for  single  degree  of  freedom  models  P  corresponds  to  both 

cree 

the  lower  and  upper  bounds  (MPCL  &  MGCL) . 

Model  C 

The  concept  of  dynamic  stability  is  similar  to  the  one  used  for  one- 


degree-of -freedom  models.  The  only  difference  is  that,  in  this  case,  there 
is  a  lower  bound  (MPCL)  and  an  upper  bound  (MGCL) .  The  lower  bound  corre¬ 
sponds  to  loads  for  which  there  is  a  possibility  of  "buckled"  motion,  while 
the  upper  bound  corresponds  to  loads  for  which  the  motion  will  definitely 
be  "buckled".  Because  the  existence  of  bounds  is  dependent  upon  the  value 
of  A,  the  discussion  will  be  based  on  the  range  of  A  values. 


Fig.  4.1.  Total  Potential  versus  Displacement  for  Various  Loads  (Model  A) 


0 


5 


10 


Fig.  A. 


radians 


2.  Static  and  Dynamic  Critical  Loads  (Model  A). 
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Fig.  4.3.  Static  and  Dynamic  Critical  Loads  (Model  B). 
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(a)  For  A  <  3  there  are  no  saddle  points  and  the  system  behaves  as  a  one- 
degree  -of  -freedom  system.  Therefore,  for  this  case  the  governing  equations 
are 

-  0  ;  (A-l-r2)  r  =  Q; 

2-P 

and  d  U_ 

- r- <  0  at  the  solution  (4.5) 

dr 

Note  that,  for  this  case,  "buckled"  motion  is  both  possible  and  guar¬ 
anteed.  Thus  the  upper  and  lower  bounds  are  coincident. 

(b)  For  A  >  3,  there  exist  "saddle"  points  on  the  total  potential,  in  addi¬ 
tion  to  the  relative  maximum  point.  All  three  are  unstable  static  equilibrium 
points  (stationary  points).  Therefore,  the  motion  can  become  "buckled"  either 
through  a  saddle  point  or  through  the  relative  maximum  point.  It  will  next 

be  shown  that,  in  certain  instances,  "buckled"  motion  is  possible  and  in  others 
guaranteed.  Before  proceeding  with  the  analysis  the  following  observations 
are  made : 

Parenthesis  1:  For  a  specified  load  the  total  potential  has  a  lower  value 

at  a  "saddle"  point  than  at  a  relative  maximum  point.  This  can  be  proven 

-P 

through  computation  of  at  the  corresponding  points. 

Parenthesis  2:  Regardless  of  the  value  of  p  and  A  (for  A  >  3),  the  total 
potential  at  any  s  position,  but  a  fixed  r-position  (r  •  r  <  Jh  -  3)  has  a 
higher  potential  than  that  of  a  position  characterized  by  r  -  r  and  an  s- 
poaltion  on  the  ellipse 

r2  +  3s2  -  A  -  3  (4.6) 

Note  that  the  ellipse,  Eq.  4.6,  defines  the  locus  of  static  equilibrium 
points  ("saddle"  points). 
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Proof:  Let  AU^,  denote  the  difference  In  total  potential,  between  the  any 

8>position  and  that  of  an  s-poaition  on  the  ellipse  (r  *  r).  Than,  from 
Eqs .  (2.11  and  4.6) 

Au£  -  [r2  +  9s2  -iJK  r  +  A  +  £  (A  -  r2  -  3s2)2  -  ZpC^A  -  r)'j 

-[4A  -  -  2r2  -  2 JK  r  -  2p(v/A  -  r)] 

-2  2  9  1  -2  2  2 

=  3  (r  +  3s  -  A)  +  J  +  J  (A  -  r  -  3s  V 

=  ~  (A  -  3  -  r2  -  3s2)2s  0  Q.  E.  D.  (4.7) 

Clearly,  the  difference  is  zero  when  s  is  on  the  ellipse  and  positive  for 
all  other  a. 

Parenthesis  3:  For  p  £  JK  +  2^t-3[loads  higher  than  the  static  critical 
load  -  see  Eq.  (2.20)],  if  r^  <  r^  i  ./A  -  3,  then  the  total  potential  on  the 
ellipse,  Eq.  (4.6),  is  higher  at  r^  than  at  r^. 

Proof: 

Let  Au£  -  U,J!(r2)  -  U^rj)  (4.8) 

Then,  by  Eqs.  (2.11)  and  (4.8) 

AU.jp  =  2 (i'2_r^)  [p  -  >/A  -  (r2  +  rt)  ]  (4.9) 

Since  r2  -  r^  is  positive  A'J.j!  i8  positive  if  p  >  +  (r2  +  r^). 

But  p  2  JK  +•  2 A/h  -  3,  thus 

if  JK  +  ijh  -  3  >  JK  +  (r2  +  r^)  then 
definitely  p  >  vrA  +  (r2  +  r^. 

Clearly,  from  the  statement  of  the  parenthesis  2,jt\  -  3  >  r^,  +  rp  w^tch 
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conclude*  the  proof. 

A*  Che  load  ia  increased  from  zero,  at  low  values  of  p,  the  zero 
potential  lines  in  the  rs-space  enclose  only  the  near  static  equilibrium 
point  and  the  motion  is  "unbuckled".  At  some  value  of  the  load  the  first 
unstable  point (s),  at  which  the  total  potential  can  become  zero,  is  (are) 
the  "saddle"  point(s)  according  to  Parenthesis  l.  At  this  load  there  is 
a  possibility  of  "buckled"  motion  through  the  saddle  point.  This  load, 
then,  is  called  MPCL.  The  governing  equations  for  finding  this  critical 
load,  as  well  as  the  corresponding  position  (sr  coordinates)  of  the  saddle 
point  are: 

UTP  -  (r2  +  9a2  -  IJKt  +  A)  +  “  (A  -  r2  -  3s2)2  -  2pC^T-r)  -  0 

(A  ~  1  -  r2  -  3a2)r  -  p  - 

A  -  3  -  r2  -  3a2  -  0  (4.10) 

Note  that  "saddle"  points  are  unstable,  thus,  there  is  no  need  for  applying 

the  stability  requirements. 

The  solution  to  Eqs.  (4.6)  yields 

(MPCL):  PcraB  *  3<VT-  1) 

l 

and  "saddle"  point  at  r  «  JK  -  j;  s  -  ±  (J\  -  £■)  2.  (4.11) 

The  solutions  to  Eq*..  (4.5)  and  (4.10)  are  plotted  on  Fig.  (4.5). 

As  far  as  the  case  of  the  upper  bound  is  concerned  (MGCL) ,  there  are 
two  ways  through  which  a  guaranteed  "buckled"  motion  can  be  achieved.  One 

•  D 

way  is  to  require  U  at  the  relative  maximum  unstable  static  point  to  be  zero 


1 


0 


In  this  case  the  motion  is  definitely  buckled  and  the  critical  load  can 
be  obtained  from  the  solution  of  Kqs .  (4.5).  The  second  way  is  for  loads 
which  are  equal  to  static  critical  load  for  asymmetric  buckling  (JK  4-  /A  -  3 ) . 
In  this  case,  although  the  total  potential  at  the  relative  maximum  point  is 
higher  than  zero,  guaranteed  "buckled"  motion  can  be  achieved  because  of 
Parentheses  2  and  3.  In  this  latter  case,  P  (MCCL)  is  given  by  the  expres¬ 
sion  for  the  static  load.  Both  results  are  shown  graphically  on  Fig.  4.5. 

Note,  that  the  upper  bound  is  the  smallest  load  computed,  either  by 
the  solutions  of  Eqs.  (4.5)  (for  s  A  <  7.3),  or  by  the  static  critical 
load  (for  A  ^  7.3).  Moreover,  for  very  large  A-values  the  upper  and  lower 
bounds  approach  each  other  (3 vJK) . 
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SECTION  V 


CONSTANT  LOAD  OF  FINITE  DURATION 

Statement  of  the  Problem 

Consider  a  model  at  its  natural  (unloaded)  position.  At  time  t  =  0 
a  constant  load  P  is  applied  suddenly  on  the  system  and  it  acts  only  for  a 
finite  duration  time  t»  T  .  After  the  release  of  the  force  P,  the  system 
moves  because  of  the  imparted  energy  during  the  action  of  the  load  P  .  The 
concept  of  dynamic  stability  for  this  particular  load  case  is  similar  to 
the  concept  used  in  the  cases  of  "ideal  impulse"  (Section  III)  and  of  "con¬ 
stant  load  of  infinite  duration"  (Section  IV).  Thus,  for  all  three  cases  the 
concept  of  dynamic  stability  is  based  on  the  definition  of  "buckled"  and 
"unbuckled"  motion. 

If  the  energy  imparted  into  the  system,  through  the  load  P,  is 

insufficient  for  the  system  to  reach  the  unstable  static  equilibrium  point 

on  the  "zero  load"  total  potential,  with  zero  velocity  (zero  kinetic  energy) 

the  motion  is  called  "unbuckled".  In  this  case,  the  system  is  dynamically 

stable.  Consequently,  the  criterion  for  dynamic  stability  requires  that 

the  dynamically  stable  system  possess  total  energy  (at  the  release  time,  T  ) 

at  a  level  below  the  level  of  the  potential  at  its  unstable  static  equilibrium 

point  for  the  "zero-load"  system.  This  is  fully  substantiated  by  mathematical 

arguments  in  the  next  section.  On  each  individual  model,  the  criterion  is 

invoked  and  estimates  for  critical  conditions  are  found.  The  extreme  cases 

of  T  -•  (constant  load  of  Infinite  duration)  and  T  -»  o  (ideal  impulse)  are 
o  o 

special  cases  of  the  present  one. 

General  Procedure 

Since  the  system  is  conservative  during  the  action  of  the  load,  then 
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(5.1) 


u!  f  TP  =  0  for  0  t  T 
T  o 


where  Tq  is  the  time  of  release  of  the  load  P.  Similarly, 


l£  +  T°  =  U°  (T  )  +  T°  (T  >  for  t  >  T 
I  1  o  o  o 


(5.2) 


The  continuity  of  the  kinetic  energy  at  the  time  of  release,  T  ,  is  expressed 

o 


T°  (T  )  =  TP  (T  ) 
o  o 


(5.3) 


Use  of  Eqs.  (5.1)  and  (5.3)  into  Eq.  (5.2)  yields 


U°  +  T°  -  U°  (T  )  -  UP  (T  )  for  t  ;r  T 
i  I  O  I  o  o 


(5.4) 


If  Lu' indicates  the  position  of  the  unstable  static  equilibrium  point  for 

P  ■  0,  the  critical  condition  is  met  (buckled  motion  Is  possible),  if  the 

load  P,  acting  for  time  To>  imparts  sufficient  energy  into  the  system  to 

equal  the  potential  U°  (L°)  of  the  "zero-load"  system  at  the  unstable  position 
o 

Lu«  Thus,  the  stability  criterion  is  expressed  by 


u°  <T0>  -  <  <TC>  *  U°  a °> 


(5.5) 


The  equality  sign  implies  a  critical  condition,  wh'le  the  inequality  refers  to 
a  dynamically  stable  situation. 

Note  that  Eq.  (5.5)  relates  the  applied  load  P  and  the  release  position, 

Lp  (one  equation  in  two  unknown  quantities,  if  one  views  the  release  position 

Lp  at  time  Tq,  as  one  quantity).  In  addition,  through  equation  (5.1),  one  may 

relate  the  applied  load,  P,  the  release  position,  L  ,  and  the  release  time,  T  , 

P  o 

for  a  specified  path  of  motion.  For  one-degree-of-frecdom  systems  there  is 
only  one  path,  but  for  two-degreo-of- freedom  systems  there  are  several  paths; 


t 


for  such  cases,  the  path  is  determined  through  the  brachistochrone  problem 

of  the  system.  This  is  so,  because  one  may  view  the  problem  as  follows: 

Find  the  smallest  time  T  for  a  given  load  P  such  that  the  motion  can  become 

o 

"buckled". 

A  critical  condition  is  obtained  by  specifying  the  load  and  finding, 
through  the  simultaneous  solution  of  the  above  two  mentioned  equations,  the 
corresponding  release  time,  T 

cr 

These  steps  are  next  clearly  demonstrated  for  each  one  of  the  three  models. 
Model  A 


For  this  particular  model,  the  unstable  stationary  points  on  the  "zero-load" 
total  potential  are  located  at  0  =  +  Tt/2  [sec  Eq.  (2.2)].  Then,  Eq.  (5.5), 
with  the  equality  sign),  becomes 


p  (cos  0  -  cos  0  )  -  [Jl  -  yi+sTn©-]2 

o  cr  ^  o 


(5.6) 


where  0  is  the  release  position,  L  ,  of  the  system.  Note  that  0  =  0(T  ). 

cr  P  cr  v  ocr/ 

This  equation,  Eq.  (5.6),  is  one  of  the  two  needed  equations  for  finding  a 
critical  condition.  Before  finding  the  second  equation,  from  Eq.  (5.1),  the 
following  simplifications  arc  introduced. 

First,  it  is  assumed  that  the  bars  of  the  system  are  weightless  and  the 
mass,  m,  of  the  system  is  concentrated  at  joint  B  (see  Fig.  2.1).  With  this, 
the  kinetic  energy  of  the  system  is  given  by 


„  1  ,2,d0,2 

T  “  2  mL  (dt) 


(5.7) 


Next,  the  following  nondimens ionalized  time  parameter,  t,  is  introduced. 


T  »  t(~)* 


(5.8) 


Finally,  the  kinetic  energy,  T,  is  nondimensionalized  in  a  manner  similar 
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to  the  total  potential 


T  =  T/kL 


(5.9) 


Then  Eq.  (5.1)  becomes 


— p  -p  -p  HO 

1C  +  T  =  lC  +  (t-) 

T  T  dT 


(5.10) 


The  critical  condition  is  obtained  from  Eqs.  (5.6)  and  (5.10).  The  approach 
one  may  use  is  as  follows:  first  specify  the  value  of  p  (these  values  must 
be  larger  than  the  critical  load  corresponding  to  the  infinite  duration  case) 
and  through  Eq.  (5.6)  solve  for  0  ;  then  use  Eq.  (5.10)  to  find  the  correspond¬ 

ing  duration  time  t  ,  or 


r  „  r'cr 

o  •'  0  -s/ItjP 
cr  n  Urr 


(5.11) 


where  is  given  by  Eq.  (2.2). 

Note  that  computationally  it  is  much  simpler  to  assign  values  of  0 

cr 

(starting  with  0q  +  0^,  where  0^  is  a  small  increment  and  increase  it  to  values 

approaching  the  angle  corresponding  to  the  infinite  duration  case)  and  solve 

for  p  and  t  from  Eqs.  (5.6)  and  (5.11)  respectively. 

°cr 

A  computer  program  is  written  and  data  arc  obtained  through  the  Georgia 
Tech  high  speed  digital  computer  CYBER  70,  Model  74-28.  These  data  are  present¬ 
ed  graphically  on  Figs.  5.1  and  5.2.  Fig.  5.1  shows  critical  condition,  for 

various  imperfection  angles,  0  ,  in  terms  of  p  and  t_  .  Note  that  as  the 

o  °cr 

critical  release  time,  t  ,  is  increased,  the  corresponding  load  approaches 

°cr 

the  value  of  pcr  Cor  the  infinite  duration  case.  The  same  data  are  presented 

on  Fig.  5.2  as  plots  of  (pr  ) _  versus  t_  .  Through  this  figure,  it  is 

o  cr  wcr 

shown  that  as  t  -*  o  the  values  of  (pr  )  approach  the  critical  impulse  values, 
*•  o  c  r 
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Impulse,  (PT0) »  versus  critical  Duration  Time,  ,  (Model  A) 


as  obtained  in  Section  Ill. 


Model  B 

For  this  particular  model  the  unstable  stationary  points  on  the  ’’zero-load" 
total  potential  are  also  located  at  6  =  +  n/2  [see  Eq.  (2.5)].  Then,  Eq.  (5.5) 
becomes 

2  p  (1  -  cos  0cr  +  e  sinOcr)  =  1  (5.12) 

_  2 

where  e  =  e/L,  p  *  PL/ka  and  is  the  position  of  the  system  at  the  critical 
time  of  release,  r 

°cr 

According  to  the  nondimensionalization  employed  for  this  model  in  Section 
III,  the  kinetic  energy  and  time  parameters  are 

T  =  02  ,  t  =  (ka2/l) *t  (5.13) 

where 

0  *  dO/dt. 


For  this  model  then,  Eq.  (5.1)  becomes 


(5.14) 


From  Eq.  (5.14)  one  may  solve  for  the  crital  duration  time,  t  ,  or 

°cr 


T 


cr 


(5.15) 


where  the  expression  for  u£  is  given  by  Eq.  (2.5).  For  this  model  as  well  as 

for  model  A,  a  computer  program  is  written  and  data  are  generated  for  finding 

critical  conditions  (p  and  corresponding  t  ).  Computationally,  values  of 

cr 

0cr  are  assigned  and  the  value  of  p  is  found  from  Eq.  (5.12).  Then,  the 

critical  release  time  rQ  is  evaluated  from  Eq.  (5.15).  The  generated  data 

cr 

are  presented  graphically  on  Figs.  5.3  and  5.4.  Note  that,  as  the  duration 


time  becomes  large,  the  corresponding  load  approaches  the  critical  value 


for  the  case  of  constant  load  of  infinite  duration,  Pcr  (see  Fig.  5.3). 

00 

Similarly,  as  the  duration  time  approaches  zero  the  corresponding  (prQ)cr 
approaches  the  value  of  the  critical  impulse  (see  Fig.  5.4).  This  value 
is  equal  to  1/e  and  it  can  also  be  derived  analytically  from  Eqs.  (5.12) 
and  (5.15).  If  the  expression  for  is  used  in  Eq.  (5.15),  one  obtains 

0 _  dO 

°cr  Jo  JlZESZjEIM  77^7  (5.16) 

1  -  cos  8  +  S  sin® 

cr  cr 

Next  it  is  observed  that  as  the  critical  duration  (release)  time 
becomes  very  small  the  corresponding  release  position,  8cr,  also  becomes 
very  small,  for  every  fixed  value  of  e.  Then,  for  very  small  8cr  values; 
and  consequently  8-values,  Eq.  (5.16)  can  be  written  as 


P9cr  d0  ,  o  P 

J-  (9,J  cr.'o  (6Cr  °cr 


cr,L_^i 


cr 


cr 


20cr  <r->* 
cr 


Cr  =  28 
o  cr 


Similarly,  from  Eq.  (5.12) 


(5.17) 


Thus, 


«Vcr 


1/e. 


(5.18) 


Model  C 

The  expression  for  the  "zero  load"  total  potential  is  given  by 

Uj  -  (r+  9s2  -  Z/Ar+A)+  %  (A-r2-3s2)2  (5.19) 

Starting  with  ,  requiring  equilibrium  and  performing  a  static  stability 


analysis  the  following  stationary  points  arc  obtained  (see  Section  III). 


Stable(Relative  min) 
Unstable  (Relative  max) 
Stable  (Relative  min) 

Unstable  (Saddle  point) 


Pt.  1  at 

Pt.  2  atj^fi  -  4  )/2,o] 

Pt.  3  a t[(JK  +  )/2,o] 

Pt.  4  at  [Va/2,  a/a-4/2J 
Pt.  5  at  [VA/2,  -  V^-4/2]  Unstable  (Saddle  point) 

For  simplicity,  it  is  assumed  that  the  three  bars  are  weightless  and 
the  only  massive  parts  of  the  system  are  two  masses,  m,  concentrated  at 
the  joints  B  and  C.  Moreover,  the  nondimensionalized  kinetic  energy  is 

—  x 

expressed  by  X  =  -75 — r  ,  where  T  characterizes  the  kinetic  energy  and 

0^kl/ 


X  =  i  ^‘2;=  [1+S'27^V 

T  2  L\dT;  +  '  dT.  I  2  L  \  dT / 


(5.20) 


where 


(0-& 


and  t 


Clearly,  saddle  points  exist  for  A  >  4.  For  this  range  of  A-values,  the 
"zero  load"  total  potential  value  at  the  saddle  points,  pts.  4  &  5,  is 
smaller  than  the  corresponding  value  at  the  relative  maximum,  pt.  2.  On 
the  basis  of  this  observation  the  motion  can  possibly  become  "buckled" 
through  the  "saddle"  points,  pts.  4  &  5.  The  corresponding  condition  for 
this  case  is  a"possible  critical  condition".  On  the  other  hand  if  the 
imparted  energy,  by  the  applied  force  at  the  release  time,  is  sufficient  to 
reach  the  relative  maximum  (unstable)  static  equilibrium  point,  pt.  2, 
"buckled"  motion  is  guaranteed  and  the  corresponding  critical  condition  is 
a  guaranteed  one.  Xhe  former  is  termed  sufficient  condition  for  dynamic 
stability,  while  the  latter  sufficient  condition  for  dynamic  instability  by 
Hsu  [173. 
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Next,  the  computational  procedure  for  finding  the  possible  critical 


condition  is  outlined. 

Through  Eq.(5.5)(in  a  nonditnensionalized  form)  one  obtains 

2p  fc/K-r)l  =  \  (A-l)  (5.21) 

T=To 


where  t  is  the  release  time, 
o 

Moreover,  Eq.  (5.1)  in  nondimensionalized  form  holds  for  (0  <  T  <  tq). 

For  a  given  path  of  motion,  integration  of  Eq.  (5.1),  yields  a  relation 

between  the  time  of  release  and  the  position  at  that  instant.  Note,  that 

the  problem  has  been  cast  in  the  following  terms:  for  a  given  load,  p, 

find  the  smallest  release  time,  t  ,  such  that  the  system  may  reach  an  - 

cr 

unstable  point  (saddle  point  for  the  minimum  possible  critical  condition) 

with  zero  velocity,  Eq.  (5.21).  Since  one  is  Interested  in  obtaining  the 

smallest  release  time,  t  ,  and  since  the  position  at  the  time  of  release 

°cr 

is  path  dependent,  one  can  solve  the  problem  by  considering  the  associated 

brachistochrone  problem.  The  brachistochrone  problem  makes  use  of  Eq.  (5.1) 

for  this  system,  and  through  its  solution  one  obtains  the  relation  between 

the  smallest  release  time,  Tocr,  the  position  at  the  instant  of  release,  as 

well  as  the  path  that  yields  T  .  The  details  of  the  solution  to  this 

°cr 

brachistochrone  problem  are  presented  in  Appendix  A.  The  solution  to  the 


brachistochrone  problem  yields  that  the  path  is  characterized  by  s^o  and  the 
relation  between  TCcr  and  the  position  of  the  system,  rcr,  at  tq  is 


A;(vft-r)-2(rVA)2-(A-r2)2 


(5.22) 


Computationally,  it  is  simpler  for  one  to  assign  values  of  rcf  (starting  with 
values  close  to  the  initial  position,  r  =  JR  and  s  =*  0),  solve  for  p  through 


Eq.  (5.21)  and  then  for  TQ  through  Eq.  (5.22). 


Note  that  for  the  case  of  the  minimum  guaranteed  critical  condition 

Eq.  (5.21)  is  replaced  by  a  comparable  equation  which  employs  the  value  of 

the 'kero  load"  total  potential  at  the  relative  maximum  unstable  static  point. 

Numerical  results  are  presented  graphically  on  Figs.  5.5  and  5.6,  for  the 

minimum  possible  critical  condition  only,  and  various  values  of  A.  The  curves 

of  Fig.  5.5  depict  critical  conditions  in  terms  of  applied  load,  p,  versus 

critical  release  time,  T0  .  One  may  observe  that  as  the  T  increases,  the 

cr  °cr 

corresponding  load  approaches,  asymptotically,  the  value  of  pcj_  for  the  infinite 

duration  time.  Fig.  5.6  presents  the  same  results  as  Fig.  5.  5,  but  in  terms 

of  (pr0)cr  versus  critical  release  time  tq  .  Note  that  as  T0  approaches 

cr  cr 

zero,  the  value  of  (pr0)cr  approaches  that  of  the  critical  ideal  impulse. 


Impulse,  <ro)„  versus  critical  Duration  Time,  t  ,  (Model  C>. 


SECTION  VI 


THE  INFLUENCE  OF  PRELOADING- CONSTANT  LOAD  OF  FINITE  DURATION 

Statement  of  the  Problem 

PQ 

Consider  a  model  at  its  stable  equilibrium  position  Lg  ,  when 

subjected  to  an  initial  static  load  Pq.  At  time  t  =  0  an  additional 

constant  load  P  is  suddenly  applied  to  the  system  and  acts  only  for  finite 

duration  time  t  *  T  .  After  the  release  of  the  force  P,  the  system  moves 

o 

because  of  the  acquired  total  energy  during  the  action  of  the  load  P. 

The  system  will  be  called  dynamically  stable  if  its  motion  is  "unbuckled", 
in  the  sense  described  in  previous  sections.  Since  the  systems  under 
consideration  exhibit  either  limit  or  unstable  bifurcation  point  insta¬ 
bility,  the  system  is  stable  if  the  energy,  imparted  through  the  action  of 
the  load  P,  is  insufficient  for  the  system  to  reach  the  unstable  static 
equilibrium  point  for  "P^load"  total  potential  of  the  system  with  zero 
velocity  (zero  kinetic  energy).  For  each  individual  model,  the  criterion 
la  invoked  and  estimates  for  critical  conditions  are  found.  The  extreme 
cases  of  To~  •  (constant  load  of  infinite  duration)  and  Tq—  0  (Ideal  Im¬ 
pulse)  are  treated  as  special  cases. 

General  Procedure 

The  concept  of  Dynamic  Stability  and  the  general  procedure  as  well  are 
extensions  of  those  used  in  the  cases  of  the  constant  load  of  finite 
duration.  The  equilibrium  positions  of; the  preloaded  system  are  given 
as  solutions  to 

P 


where  L  is  the  position  of  the  system.  Thus,  one  may  find  all  the  "P  - 

o 

load"  static  equilibrium  positions  including  the  near  stable  position 
P  Pfl 

L  as  well  as  the  unstable  position  L  ,  through  which  "buckled"  motion 

can  be  realized  (sec  Fig.  6.1). 

Keeping  the  same  generalized  coordinates  for  all  models  and  the 
same  expressions  for  the  total  potential  and  for  the  kinetic  energy,  one 
may  apply  the  concepts  already  developed.  These  are  next  explained  through 
the  use  of  Fig.  6.1,  which  holds  for  onc-degree-of- freedom  systems,  but 
the  explanation  is  applicable  to  all  f inite-degree-of-freedun  systems. 


The  system  is  initially  loaded  quas i-s tatically  by  load  p  and  it 
Pq  ° 

reaches  point  A  (L  *  ;  stable  static  equilibrium  point).  Then,  a  load 

p  is  applied  suddenly  for  a  time  tq  (finite  duration).  At  the  load  p 

is  removed. 

P 

-P  _  p  -  O 

A  potential  U  is  defined,  such  that  U  =  U  at  t  =  0  or 

P  *  *  A 

L  *  Lg°  (see  Fig.  6.1). 


P  +P  ,  P  . 


n  *  _  '  *  r  »  /  *  \  A  •  a  /  A  \ 

0T  '  DT  +  [0T  («*#)  -  °T°  (LS°)J 


(6.2) 


Since  the  system  is  conservative,  then  during  the  action  of  p  one  may 

write 

0?  ♦  iP°+P  -  £  (£)  !  0  *  T  •  T„ 

_p  +p 

where  T  is  the  kinetic  energy  of  the  system.  This  equation,  Eq.  (6.3), 
is  equivalent  to  Eq.  (5.1)  of  Section  V.  Making  use  of  Eq.  (6.2),  Eq.  (6.3) 


becomes 


P  +P  P  +P  P  +P  .  P  . 

0T°  +  *  °  '  ®T°  (lS°)  !  °  ’  '■  ’o 


(6.4) 
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For  times  greater  than  t  ,  the  system  is  also  conservative  and  conser¬ 
vation  of  energy  yields 


V  V  1* 

UP  +  T  °  =  UT°  <to)  +  T  °  (tq)  ;  T  >  tq 


(6.5) 


where  the  right-hand-side  represents  the  level  of  the  system  energy  at  the 

release  time  t  ,  and  Eq.  (6.5)  is  equivalent  to  Eq.  (5.2)  of  Section  V. 

If  the  force  p  has  imparted  sufficient  energy  into  the  system,  such 

P 

that  it  can  reach  the  unstable  point  B  (L  =  L^0;  see  Fig.  6.1)  on  the 
"p  -load"  potential  with  zero  kinetic  energy  (velocity),  then  "buckled" 
motion  is  possible,  and  the  system  becomes  dynamically  unstable. 

The  governing  equations  for  predicting  critical  conditions  are  obtained 
from  Eqs.  (6.4)  and  (6.5). 


By  requiring  kinematic  continuity  at  i  =  tq  one  may  write 

P  +P  P 

T  °  (to)  =  T  °  (to) 


(6.6) 


From  Eq.  (6.4)  one  may  write 


P  +P  P  +P  ,  P  „  P  +P 
Z.  o  _  -  o  f,o  -  o 


/  »  \  1  'l 

•  V  <To> 


(6.7) 


Substitution  of  Eqs.  (6.6)  and  (6.7)  into  Eq.  (6.5)  yields 


UT°  +  T 


*  P  P  +P  P  +P,  P  \ 

°  -«x°  ('o>  -  GX°  (To>  +0T°  (LS/ 


(6.8) 


for  r  >  t 


Next,  a  critical  condition  exists,  if  the  system  can  reach  point  B  (see 

r  Po  (  Po^  n 

Fig.  6.1)  with  zero  kinetic  energy  T  °  j  *  OJ.  Thus,  denoting  the 


critical  condition  by  p,  t  ,  one  may  write 

cr 


P  +P,  P 


*  /  \  r  ^1/  \  a  /  *  \ 

U_°  (l  °  !  =  U_°  (t  j  -  V  °  (r  I  +  0°  (lc°) 

T\u/  T\o  /  T  \  o  '  T  \  S  > 


(6.9) 
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p 

ts 

V 

l! 


Note  that  Eq.  (6.9)  relates  p,  Tq  and  the  position  of  the  system, 

cr 

L  ,  at  the  Instant  of  the  release  of  the  force  p.  Please  observe  (as 
cr 

also  explained  in  Section  V)  that  the  above  described  critical  condition 
depends  on  the  evaluation  of  two  parameters,  p  and  t  .  One  approach  is 
to  prescribe  and  find  the  corresponding  pc^  and  the  other  is  to  pre¬ 
scribe  p  and  find  the  corresponding  t  .  The  two  are  equivalent.  Re- 

cr 

gardless  of  the  approach,  Eq .  (6.9)  relates  three  parameters  p,  r  (or 

cr 

p  ,  Tq)  and  L  (the  position  of  the  system  taken  as  one  parameter). 

The  second  (needed)  equation  is  obtained  from  Eq.  (6.4).  This  equa¬ 
tion  (see  Section  V)  is  used  to  relate  the  load  p,  the  time  of  release, 

t  ,  and  the  position  of  the  system  a*  the  instant  of  release.  In  order 
o 

to  find,  for  a  prescribed  load,  the  position  of  the  system  at  the  instant 
of  release,  one  must  specify  the  path  of  motion.  For  one-degree-of- freedom 
systems  there  is  only  one  path  and  this  is  easily  accomplished  [see  Eqs . 
(5.13)  and  (5.17)].  On  the  other  hand,  for  a  mu lti-degree-of- freedom 
system  there  are  numerous  possible  paths  leading  to  a  multitude  of  posi¬ 
tions  for  a  given  release  time.  In  such  cases,  if  one  is  interested  in 
finding  a  lower  bound  for  the  critical  conditions,  he  may  find  the  path 
that  yields  the  smallest  possible  time.  Thi°  may  be  accomplished  by 
solving  the  corresponding  "brachistochrone"  problem  (see  Appendix  A).  Hie 
solution  to  the  "brachistochrone"  problem  yields  the  path  of  motion  which 
yields  the  smallest  possible  time.  Thus,  Eq.  (6.4)  along  with  the  path 
of  motion,  relates  the  release  time  tq  and  system  position  at  the  instant 
of  release  [see  Eq.  (5.24)].  These  steps  are  clearly  demonstrated  for 
each  of  the  three  models,  in  the  subsequent  articles. 

Moreover,  if  one  is  dealing  with  a  deflection  limited  design  (the 
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1*0  -P  P 

position  of  tho  system  cannot  exceed  T.  ‘  -  see  Fig.  6.1),  then  U  °(I.  °) 

1.  T  u 

-F  P 

is  replaced  by  U  o(L  o) ,  in  the  outlined  computational  procedure. 

1  la 

Parenthesis :  The  cases  of  ideal  impulse  and  suddenly  applied  load  of 
constant  magnitude  and  infinite  duration,  may  be  obtained  as  special  cases 
qtf  the  present  procedure. 

However,  critical  conditions  for  these  two  load  cases  may  also  be  ob¬ 
tained  independently. 

For  the  ideal  impulse  case,  one  may  relate  the  impulse  to  an  initial 
kinetic  energy,  and  from  conservation  of  energy 


GPo  +  *Po  „  qPo  (1P0>  + 

T  T  's'  l 


(6.10) 


Then  TP°  is  critical  (related  to  Lhe  critical  ideal  impulse)  if  the 
system  reaches  position  LP°  (pt  B  of  Fig.  6.1)  with  zero  kinetic  energy. 
Thus 


TP°  =  UP°(l.Po)  -  UPo  (lf°) 
l  T  s  T  u 


(6.11) 


For  the  second  extreme  case  (T0*")>  P  may  be  obtained  from  Eq.  (6.3) 
or  (6.4),  which,  for  this  ease,  holds  true  for  all  t  (0  v  t  -■  °°).  Thus,  P 

c 

corresponds  to  the  solution  of 


P  +  PP+P  P+P 

0  °  (L  °u  )  =  U  °  (l.Po)  (6.12) 

Po  +  p 

where  L  denotes  the  unstable  static  equilibrium  position  of  the 

system  when  the  static  load  is  equal  to  +  p. 
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Model  A 


For  this  particular  model,  the  static  stability  analysis  is  presented 
in  Section  V.  The  geometry  of  the  model  is  given  in  Fig.  2.1  and  the  static 
response  in  Fig.  2.2. 

In  evaluating  the  effect  of  preloading  three  imperfection  angles  are 
chosen  (0q  *  0.005,  0.010,  0.020),  and  for  each  0q  -  value  the  system  is 

initially  loaded  quasis tatically  with  a  pQ-load  smaller  than  the  p^-static. 
Then,  the  system  is  loaded  dynamically.  The  following  values  are  used  in 
the  dynamic  analysis: 


e 

=  0.005 

;  p 

=  0.440 

9 

P 

*  0.340,  0.380, 

and 

0.420 

o 

cr 

o 

e 

-  0.010 

;  p 

*  0.415 

• 

9 

P„ 

-  0.300,  0.360, 

and 

0.400 

o 

rcr 

ro 

0 

-  0.020 

:  pcr 

-  0.384 

9 

po 

*  0.300,  0.340, 

and 

0.370 

First,  the  extreme  cases  (Tq  -*  o  and  -*  <»)  are  analyzed,  by  employing 
Eqs  (6.11)  and  (6.12)  respectively. 

The  ideal  impulse,  (PrrQ)  is  related  to  the  initial  kinetic  energy 
tin  the  nondimens ionalized  form;  see  Eqs.  (3.8)  and  (3.9)]  by  the  expression 


(pO  = 


sin(9 


sft'-f 


(6.13) 


where  9^  is  the  stable  static  position  (angle  9)  under  pQ  load. 

The  critical  ideal  impulse  (pt  )  ,  is  obtained  by  substituting  Eq. 

o  cr 

(6.11)  into  Eq.  (6.13),  or 


OVcr  '  71V, 

“T.  ) 


[<%»:)  ■ 


-P„,  P,, 


(6.14) 
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where  0  °  is  the  unstable  static  position  under  load,  and  the  expression 
for  the  total  potential  is  given  by  Eq.  (2.2)  by  using  pQ  wherever  p  appears 
in  the  equation.  The  numerical  results  for  all  0o>  pq  combinations  are 
presented  in  tabular  form  on  Table  6.1. 


Table  6.1  Critical  Ideal  Impulse,  (PT0)cr>  Model  A. 


0  * 
o 

3.005 

0  =  ( 
o 

3.010 

CD 

o 

N 

O 

020 

o 

Q. 

(PT  ) 

'r  o  cr 

% 

''"o'er 

o 

a. 

(pr  )rr 
o  cr 

0 

162.000 

0 

81.000 

0 

40.000 

0.340 

8.246 

0.300 

6.980 

0.300 

2.900 

0.380 

3.498 

0.360 

2.247 

0.340 

1.207 

0.420 

0.758 

0.400 

0.486 

0.370 

0.338 

0.440 

0 

0.416 

0 

0.384 

0 

Note  that  the  first  row  results  are  obtained  t  rom  Pq.  0.10)  of 
Section  III.  Note  also  that  as  the  value  ol  p  approaches  the  value  of 
the  static  critical  load  the  additionally  imposed  (critical)  impulse  tends 
to  zero.  This  is  reflected  hy  the  results  ot  the  last  row  (Table  6.1). 

The  critical  load  for  the  case  ol  r  -  »,  p  ,  Is  obtained  the 

following  steps,  for  a  given  0^,  p^  combination  (a)  solve  Eq.  (2.3)  given 

p  1  ' 

below,  for  0^  (stable  position) 

PQ  *  [./l  +  sin  8  -  ./l  VTinT*  j  cot  8/  ./ 1  V  s in  0  (6.15) 

(b)  Static  equilibrium  positions  arc  character  i/.od  by  Eq .  (2.3)  for 

loading  p  +  p,  or 
o 

p  +  py  -  [y,  +  sin”  -  J \  +  si n"lT~j  cot  0/  ./TV  sin  0  (6.16) 
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(c)  Eq.  (6.12)  for  this  model  is  given  below. 


L/l  +  sin  0  -  */l  +  sin  G  -  (p  +p)  (cos  0  -  cos  0)  = 
L  OJ  o  o 


C  \/l  +  sin  0P  -Vl+sinO-l-p  ( 

s  o  o  \ 


cos  9  -  cos  0 
o  s 


(6.17) 


P  +P 
o 


The  simultaneous  solution  of  Eqs.  (6.16)  and  (6.17)  yields  9  and 

Pcr  .  Note  that  since  Eq.  (6.17)  must  be  satisfied  at  an  unstable  point 
00 

(6Po+p)  the  stability  test  may  be  used  to  ensure  this,  or 


dV°+p 

_ T_ 

de2 


Po+P 


<  0 


I  U 

The  numerical  results  for  all  G^,  p^  combinations  are  presented  in 
tabular  form  or  Table  6.2. 

Table  6.2  -  Critical  Dynamic  Load,  p 

CO 

(Constant  Load  of  Infinite  Duration:  Model  A) 


- e 

05 - 

■■HI 

=»  U.U1U 
o 

- 9 — 

o 

-  0.020 

Po 

pcr 

OO 

P  +P 
o  rcr 

CO 

19 

Pcr 

OO 

P+P„,. 
p  cr 

OO 

■a 

pcr 

OO 

PofPcr 

OO 

0 

0.4320 

0.4320 

0 

0.4010 

0.4010 

0 

0.3690 

0.3690 

0.340 

0.0940 

0.4340 

0.300 

0.1090 

0.4090 

0.300 

0.0751 

0.3751 

0.380 

0.0560 

0.4360 

0.360 

0.0517 

0.4117 

0.340 

0.0396 

0.3796 

0.420 

0.0165 

0.4365 

0.400 

0.0136 

0.4136 

0.370 

0.0116 

0.3816 

0.440 

0 

0.4400 

0.4160 

0 

0.4160 

0.384 

0 

0.3840 

Note  that  the  first  row  results  of  Table  6.2  are  taken  from  Section  IV. 

The  results  of  the  last  row  reflect  the  fact  that  if  the  system  is  loaded 
quasistatically  up  to  the  limit  point,  then  the  additional  suddenly  applied 
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load  that  the  system  can  withstand  tends  to  zero. 

Finally,  for  the  case  of  constant  load,  p,  applied  suddenly  for  a 
finite  duration,  t  ,  critical  conditions  are  obtained  from  the  following 


steps: 


Po  Po 

(a)  From  the  static  stability  analysis  obtain  0  and  8  for  each 

s  u 


(b)  Use  of  Eq.  (6.9)  yields 


r0  — - _ _ _  _  . 

p  ^cos  Bs  -  cos  ®cr)  *  |V  1  +  sin  0^°-  "J  l  +  sin  8  J 


-  [/ 1  +  sin  0‘>o  -  V  l  +  sin  0  ] 


/  PQ  p  . 

+  p  (  cos  0  -  cos  0  °) 

o  V  u  s  / 


(6.18) 


where  6  is  the  position  0  at  the  instant  of  the  release  of  the  force  p 
cr 

(T  *  t  ). 

o 

Po 

In  Eq.  (6.1$for  a  given,  p  ,  everything  is  known  (p  ,  0  ,  0  ,  and 

p  o  o  o  s 

8  °)  except  p  and  0  .  Therefore,  Eq.  (6. 18) relates  p  and  0  for  a  critical 

u  cr  cr 

condition  to  exist. 

(c)  Since  T  °=  (^j~)  ,  then  from  Eq.(6.4)we  may  write 


d6  JV",  \  _po+l> 

d7  =  UT  (08/-UT  <0)  °r 


p  +p 

r-  0  /  P 


*  -  [V  (•  •)-  5,  (•)] 


(6.19) 


Integration  from  t  =  o  to  t  =  tq  and  use  of  the  expression  for  the  total 


potential  [see  Eq.(2. 2)]  yields 
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Note  that  Eq.  (6.20)  also  relates  6  to  p. 

cr 

A  critical  condition  is  characterized  by  (p,  t  )  that  satifies  both 
equation,  Eq.  (6. 18)  and  (6.20).  This  means  that  for  a  given  release  time, 

Tq,  find  pcr  or  for  a  given  p  find  tq  .  Computationally,  though  it  is 
easier  to  assign  values  of  9^,  solve  for  p  from  Eq.  (6. 18)  and  then  for  the 
corresponding  tq  from  Eq.  (6.20). 

A  computer  program  has  been  written  for  these  computations.  Values  of 

p 

9cr  are  assigned,  starting  with  8S°+  69,  where  69  is  very  small,  and  computing 

the  corresponding  values  of  p  and  tq  for  each  60 . 

The  results  are  presented  graphically  on  Figs.  6.2  -  6.7  for  the  three  values 

values  of  8q,  On  the  first  three  figures  critical  conditions  appear  as  plots 

of  p  versus  duration  time,  t  .  Note  that  as  T  becomes  larger  and  larger, 

°cr  °cr 

the  corresponding  value  of  p  approaches  p  (see  Table  6.2).  On  the  last 

three  figures  (6. 5-6. 7),  critical  conditions  appear  as  plots  of  (PTQ)  versus 

cr 

duration  time,  t  .  On  these  figures,  as  t  -•  o,  the  corresponding  value 
°cr  °cr 

of  (pt  )  approaches  the  critical  ideal  impulse  (see  Table  6.1). 

°  cr 
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Fig.  6.5.  Superimposed  Impulse, (pt),  versus  critical  Duration  Time,  t  ,  Preloaded  Model  A 


uperimposed  Impulse,  (pT  ),  versus  critical  Time  Duration,  t  ,  '  Preloaded  Model  A 


Model  B 


For  this  particular  model,  t lie  static,  stability  analysis  is  presented 
in  Section  II.  The  geometry  of  the  model  is  given  in  Fig.  2.3. 

In  evaluating  the  effect  of  preloading,  three  eccentricities  are  chosen 
(e  ?!  0.005,  0.010,10.020),  and  for  each  e-value  the  system  is  initially 
loaded  quasi-statically  with  a  pQ-load  smaller  than  p^-static.  Then,  the 
system  is  loaded  dynamically.  The  following  values  are  used  in  the  dynamic 
analysis: 


e 

-  0.005 

9 

Pcr 

-  0.955 

Po  =  0.35 

0.40, 

and 

0.50. 

e 

-  0.010 

9 

Pcr 

-  0.932 

Po  =  0.30 

0.40 

and 

o 

• 

o 

e 

-  0.020 

9 

Pcr 

-  0.898 

p  =  0.30 

O 

0.40, 

and 

© 

U1 

o 

First,  the  extreme  cases  (tq-h)  and  are  analyzed,  by  employing  Eqs . 

(6.11)  and. (6. 12)  respectively. 


The  ideal  impulse.  (ptq)  is  related  to  the  initial  kinetic  energy  by 
the  expression 


(PTo) 


sin 


_1 _ 

4-  e  cos 


i  /«  i  2 

l'i  J 


(6.21) 


where  6  °  is  the  stable  sialic  position  (angle  0)  under  load, 
s  ° 

The  critical  ideal  impulse,  (PT0)cr>  is  obtained  by  substituting  Eq. 
(6.11)  into  Sq •  (5.21),  or 
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(pT  ) 

r  o  cr 


r.Po/  Po\  —  /  pov-i 
Ll,T  x®.  )~  Ut(°u  ;]  <6-22> 


P  P 

sin  6S  +  e  cos  0S 


where  0^  is  the  unstable  static  position  under  pq  load,  and  the  expression 
for  the  total  potential  is  given  by  Kq.  2.5  hv  using  wherever  p  appears  in 
the  equation.  The  numerical  results  for  all  e,  p^  combinations  are  presented 
in  tabular  form  on  Table  6.3. 


Table  6.3.  Critical  ideal  Impulse,  ^PT0^{r 


e  =0.005 


(Model  (0 
e  =0.01 


e  “  0.02 


po 

<pT»’cr 

Po 

(PT  ) 
o  cr 

Po 

(PT 

0 

200 

0 

100 

0 

50 

0.35 

192 

0.30 

88 

0.30 

42 

0.40 

164 

0.40 

76 

0.40 

36 

0.50 

139 

0.50 

62 

0.50 

29 

0.955 

0 

0.932 

0 

0.898 

0 

Note  that  the  first  row  results  are  obtained  from  F.q.  (3.24)  of 
Section  III.  Note  also  that  as  (lie  value  o I  p  approaches  the  value  of 


the  static  critical  load  the  additionally  imposed  (critical)  impulse 
tends  to  zero.  This  is  reflected  by  the  results  of  the  last  row  (Table 


6.3). 


The  critical  load  for  the  case  of  t  p  »  is  obtained  by  the 

CO 

following  steps,  for  a  given  9  ,  p  combine! ion. 

P0 

a)  Solve  Eq.  (2.6)  given  below,  for  0  (stable  position) 

s 


(6.23) 


(b)  Static  equilibrium  positions  are  character ized  by  Eq.  (2.6)  for 
loading  p  +p,  or 


P  +  P 


sin  0 

o  tan  0  +  c 


(6.24) 


(c)  Eq.  (6.12)  for  this  model  is  given  below 


sin  0  -  2(pq+  p)  ( 1  -  cos  0  +  e  sin  9)  ■ 

?  P_  P_ 

sin  0  *  2  p  (l  -  cosO  +  5  sin  0  °  (6.25) 

so  s  s 

P  +  p 

The  simultaneous  solution  of  Eqs.  (6.24)  and  (6.25)  yields  0  and 

u 

p  .  Note  that  the  same  Eq.  (6.25)  must  be  satisfied  at  an  unstable  point 

OO 

PG+P 

(0u  ).  The  stability  test  may  be  used  to  ensure  this,  or 


2  Po+  P 

d  U.r° 


p  +p  <0 

.  o 
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r« 


r 


The  numerical  results  for  all  e,  p^  combinations  are  presented 
in  tabular  form  on  Table  (>.4. 

Table  6.4  -  Critical  Dynamic  Load,  p 


(Constant  Load  of  Infinite  Duration:  Model.  U) 


e  *  0.( 

1005 

e  =  0.010 

i  «  0.020 

m 

El 

P  +P 
o  *cr 

GO 

m 

pcr 

GO 

p  +p 
o  cr 

CO 

El 

Pcr 

00 

P  +P„ 
o  cr 

00 

0 

0.948 

0.948 

0 

0  912 

0.912 

0 

0.85 

0.85 

0.35 

0.60 

0.95 

0.30 

0.62 

0.92 

0.59 

0.89 

0.40 

0.54 

0.94 

0.40 

0.52 

0.92 

0.48 

0.88 

0.50 

0.45 

0.95 

0.50 

0.42 

0.92 

0.38 

0.88 

0.955 

0 

0.955 

0.932 

0 

0.932 

0.898 

0. 

0.898 

Note  that  the  first  row  results  of  Table  6.4  are  taken  from  Section  IV. 

The  results  of  the  last  row  reflect  the  fact  that  if  the  system  is  loaded 

quasi-statically  up  to  the  limit  point,  the  additional  suddenly  applied  load 

that  the  system  can  withstand  tends  to  zero. 

Finally,  for  the  case  of  constant  load,  p.  applied  suddenly  for  a 

finite  duration,  t  ,  critical  conditions  are  obtained  from  the  following  steps: 

P  P 

a)  From  the  static  stability  analysis  obtain  9  and  9  for  each  p  . 

s  u  o 

b)  Use  of  equation  (6.9)  yields 

po  .  PQ 

2p  (cos  6  -  e  sin  9  -  cos  9  4-  e  sin  9  )  ■ 

$  g  c r  cr 


.2  A, 


2  A 


sin  9  -  sin  9  •  2p  (cos  9  -  e  sin  9  c 

u  sos  s 


-  cos  9  +  e  s 


in  9P°) 

u  / 


(6.26) 


where  0  is  Ihe  post' ion  0  a*  the  ins'n  it  of  the  release  of  t  lie 
cr 

force  p  (t  »  t  ) 

In  Eq.  6.26.  for  a  given  geometry  v  and  static  load  p^,  everything 
P  p 

is  known  (p  .  e,0  °  and  0  °)  extepi  p  and  0  .  Therefore,  Eq.  6.26 

os  u  t r 

relates  p  and  8  for  a  critical  condition  to  exist. 
c  2 

c)  Since  xP+P°"(^^)  .  then  from  Eq.  6.4  one  may  write 

r.P  +P/  P«\  .Prt*P  Y2' 


dr  L  T  Vs  /  T  I 


P  +P  ,  P 


(6.27) 


Integration  from  t  =  0  to  t  =  tq  and  use  of  the  expression  for  the 
total  potential  [see  Eq.  2.5]  yields 


t  *  f  f  sin^  0  °  *  sin^  G 
o  J  *•  s 


p  p  -  T 

■2  (p  +p)  »(e  sin  9  °  -  cos  9  °  -  e  sin  0  +  cos  0)}  d0  (6.28) 

o  s  s 


Note  that  Eq.  (6.28)  also  relates  0  to  p. 

cr 

A  critical  condition  is  character ized  by  (p,  tq>  that  satisfies  both 

equations,  Eqs.  (6.26)  and  (6.28).  This  means  that  for  a  given  release 

time,  t  ,  find  p  or  for  a  given  p  find  t  .  Computationally,  though, 
o  cr  ccr 

it  is  easier  to  assign  values  of  0  ,  solve  for  p  from  Eq.  (6.26)  and  then 

for  the  corresponding  tq  from  Eq.  (6.28). 

A  computer  program  has  been  written  for  these  computations.  Values 


o 

of  0  are  assigned,  starting  with  9  +  60  •  where  6  0  ,  is  very  small,  and 

cr  8 

computing  the  corresponding  values  of  p  and  tq  for  each  60. 

The  results  are  presented  graphically  on  I'igs.  6.8  -  6.13  for  the 

three  values  of  0q.  '  On  the  first  three  figures  critical  conditions  appear 

as  plots  of  p  versus  duration  time,  t  .  Note  that  as  t  becomes  larger 

cr  °cr 

and  larger,  the  corresponding  value  of  p  approaches  p  (see  Table  6.4). 

CD 

On  the  last  three  figures  (6.11  -  6.13)  critical  conditions  appear  as  plots 

of  (dt  )  versus  duration  time,  t  •  On  these  figures,  as  t  -•  o,  the 
r  o  cr  °  o 

^  cr  cr 

corresponding  value  of  (PT0)cr  approaches  the  critical  ideal  Impulse  (see 
Table  6.3). 
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Fig.  6.10. 


Superimposed  load,  p,  versus  critical  Duration  Time,  f  ,  Preloaded 
Model  H,  e  =  0.02.  °cr 


r imposed  Impulse,  pr ,  versus  cri 


Model  C. 

For  this  particular  model,  the  static  stability  analysis  is  presented 
in  Section  II.  The  geometry  of  the  model  is  given  on  Fig.  2.4  and  the 
static  response  on  Fig.  2.6. 

In  evaluating  the  effect  of  prestress,  three  A-characteristic  coef¬ 
ficients  are  chosen  (A  =  5.0,  6.0,  8.0),  and  for  each  A“value  the  system 
is  initially  loaded  quasi-statically  with  a  pQ-load  smaller  than  the  p  - 
static.  Then,  the  system  is  loaded  dynamically.  'Jhe  following  values 
are  used  in  the  dynamic  analysis. 


A  ■ 

5.0 

9 

pcr 

-  5.154 

9 

Po 

-  2.8, 

3.4, 

4.0 

A  * 

6.0 

9 

pcr 

-  6.000 

9 

Po 

-  3.2, 

3.8, 

4.2 

A  - 

8.0 

9 

pcr 

-  7.310 

9 

Po 

*  4.8, 

5.4, 

6.0 

First,  the  extreme  cases  (tq-*  o  and  tq  -  «)  are  analyzed,  by  employ¬ 
ing  Eqs.  (6.11)  and  (6.12),  respectively. 

The  ideal  impulse.  (PTQ)*  1*  related  to  the  Initial  kinetic  energy  in 

the  nondimensionalized  form  [see  Eqs.  (3.22)]  by  the  expression 

1 

P  * 

(PTo)  -  [2  Tj°  ]2  (6.29) 


*  Q 

where  6  is  the  stable  static  position  (angle  9)  under  p  load.  The 

8  O 

critical  ideal  impulse,  (PTQ)cr>  obtained  by  substituting  Eq.  (6.11) 


into  Eq.  (6.13) ,  or 


<»Vcr  ’  *  [°t°  (*.’) 


(6.30) 


where  0y  is  Che  unstable  static  position  under  pQ  load,  and  the  expression 
for  the  total  potential  is  given  by  Eq.  (2.2)  by  using  p^  wherever  p  appears 
in  the  equation.  The  numerical  results  for  all  p^  combinations  are  pre¬ 
sented  in  tabular  form  on  Table  6.5. 


Table  6.5  -  Critical  Ideal  Impulse,  (pTu)cr > (Model  C) 


A-5.0 

A*6.0 

A=8.0 

P 

o 

IH8M 

P 

0 

warn 

P 

o 

0 

6.00 

0 

6.71 

0 

7.93 

2.8 

3.06 

3.2 

3.39 

4.8 

3.04 

mm 

2.32 

3.8 

2.67 

mm 

2.31 

HQi 

1.32 

4.2 

1.67 

6.0 

1.44 

5.154 

0 

6.0 

0 

7.31 

0 

Note  that  the  first  row  results  are  obtained  from  Eq.  (3.26)  of 
Section  III.  Note  also  that  as  the  value  of  pQ  approaches  the  value  of 
the  static  critical  load  the  addtionally  imposed  (critical)  Impulse  tends 
to  zero.  This  is  reflected  by  the  results  of  the  last  row  (Table  6.5). 


1 


The  critical  load  for  the  ease  of  T0  -  «»,  p  ,  is  oli rained  h  the 

following  steps ,  for  a  given  \,  |>u  combination 

|> 

a)  Solve  Eq.  (2.17)  given  below  lor  r^1'  (stable*  position)  and  s  =  o, 

2 


for 


PQ  K  (A-l-r^o  )  r^'  +  JR  (6.31) 

b)  Static  Saddle  equilibrium  positions  are  characterized  by  Eq.  (2.19) 
loading  PQ  +  P.  or 

p  +  Po  - 


2 


2  _  2 
r  +  3s 


A -3 


(6.32) 


c)  Eq.  (6.12)  for  this  model  is  given  below 


(r2+  9s2-  zfR.  r  +  A)  ^  j  (A-r2-3s2'2-  2(po+p)  (./A-r)- 

p  2  i  P  2  2  r-  P 

+  --(A-r  o  )  -  2P0(./A-r  «) 
s  t  s  s 


(6.33) 


Po+p 

The  simultaneous  solution  of  Eqs.  (6.32)  and  (6.33)  yields  r^  and  p  . 

W 

The  numerical  results  for  all  A,  p^  combinations  arc  presented  in  Tabular 
form  on  Table  6.6. 


Table  6.6.  Critical  Loads  (Constant  Magnitude,  Infinite  Duration  -  Model  C). 


A  «  5.0 

A  =  6. 

0 

_ 

A  -  8.( 

) 

p 

P 

OH mm 

P  - 

l'  +  I’ 

mrm 

EH 

P„  4p„,. 

ro 

cr® 

ms 

Hi 

cr„ 

o  cr« 

mm 

IQ 

°  cxia 

0 

3.70 

3.70 

0 

isa 

4.35 

■1 

5.70 

HBB1 

2.8 

1.48 

4.28 

3.2 

US 

4.88 

1.59 

1 

3.4 

1.17 

4.57 

3.8 

1.22 

3.02 

■ 

1.21 

•4.0 

.  .58 

4 . 38 

4.2 

0.85 

5.03 

HI 

0.68 

B9 

I  5.154 
1 _ 

in 

<r 

• 

i 

_ 1 

6.000 

0.0 

L_. 

JBBm 

/.  ji 

Q| 

7.  H 
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Note  that  the  first  row  results  of  Table  6.6  are  taken  from  Section  IV. 


The  results  of  the  last  row  reflect  the  fact  that  if  the  system  is  loaded! 
quasi-statically  up  to  the  limit  point,  then  the  additional  suddenly  applied 
load  that  the  system  can  withstand  tends  to  zero. 

Finally,  for  the  case  of  constant  load,  p,  applied  suddenly  for  a  finite 

duration,  t  ,  critical  conditions  are  obtained  from  the  following  steps: 

°  Po  Po  Po 

(a)  From  the  static  stability  analysis  obtain  r  and  r  s  for 

3  u  u 


each  p 


(b)  Use  of  Eq.  (6. 9)  yields 


2  2 

P  P  P  P 

2p  (r  °-  r  )  =  (r  °+  9s  0  -  z/K  r  °  +  A)  + 
s  cr  u  u  u 


,  P„P  „  ,  P  z  2  PP 

+  2  -  <r  °-  ^  >  +  2P„<  r„°- 


(6.34) 


where  r  is  the  position  r  at  the  instant  of  the  release  of  the  force  p  (t=t  ). 
cr  o 

In  Eq.  (6.34),  for  a  given  geometry.  A,  and  Static  load,  p  ,  everything  is  known 

P0  P  p0  ° 

(p  ,  e,  r  ,  r  and  S  )  except  p  and  r  .  Therefore,  Eq.  (6.34)  relates  p  and 
o  s  u  u  cr 

rcr  for  a  critical  condition  to  exist. 

-I>fpo  1  ’2  dr  2 

<c)  Since  T  *  ~(l+s  )  ( ,  then  from  Eq.  6.4  we  may  write 

P+P  r,  P_ .  P  i. 

.it-  —  °  /  -  °+  /  \  ' 


dT  •  LUI  (rs°)-,,T  (r-Aj’.  or 

.  Jo**,  P,  Jo*r, 

*-Lut  (r,  -ux 


<6.35) 


Invoking  the  same  techniques  wc  used  for  the  same  problem  but  without  prestress 
in  Section  V,  the  critical  time  tq  is  computed  on  the  symmetric  path  so. 

Integration  from  t  =  o  to  t  =  tr  and  use  of  the  expression  for  the  total 
potential  [see  Eq.  2.13]  yields 

V  pr[<r9°'  vA>2+J<A  -  rPo)2-  (r  -  -A)2-  I  (A-r)2- 


(6.36) 


Prt 

-  2  (p  +PQ)(r-  r  s°)  J  •'  dr 


Note  that  Eq.  (6.34)  also  relates  rcr  to  p. 

A  critical  condition  is  characterized  by  (p,  t  )  that  satifies  both 

equation,  Eqs.  (6.34)  and  (6.36).  This  means  that  for  a  given  release  time, 

t  ,  find  p  or  for  a  given  p  find  t  .  Computationally,  though  it  is  easier 

cr 

to  assign  values  °^rcr»  solve  for  p  from  Eq.  (6.34)  and  then  for  the  correr 
sponding  tq  from  Eq.  (6.36). 

A  computer  program  has  been  written  for  these  computations.  Values  of 

p 

r  are  assigned,  starting  with  rs°+  ®r»  where <fr  is  very  small,  and  computing 

the  corresponding  values  of  p  and  for  each  <fr. 

The  results  are  presented  graphically  on  Figs.  6.14  -  5.17  for  the  three 

values  of  A.  On  the  first  three  figures,  critical  conditions  appear  as  plots  of 

p  versus  duration  time,  t  .  Note  that  as  t  becomes  larger  and  larger,  the 

°cr  °cr 

corresponding  value  of  p  approaches  pcr  (see  Table  6.6).  On  the  last  three 

CO 

figures  (6.17  -  6.79)  critical  conditions  appear  as  plots  of  (prQ)  f  versus 

duration  time,  T  .  On  these  figures,  as  t  -*  o,  the  corresponding  value  of 
°cr  cr 

(pT  )  approaches  the  critical  ideal  impulse  (see  Table  6.5). 
cr 
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Fig.  6.14.  constant  load,  p,  versus 


Fig.  6.15.  Constant  load,  p,  versus  Critical  Duration  Time,  tq  ,  preloaded  Model  C,  A  =  6.0 


tnpulse,  (pT  ),  versus  critical  Dur;  :ion  Time,  ~q  ,  preloaded  Model  C 
°  cr 
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SECTION  VII 


THE  INFLUENCE  OF  SMALL  DAMPING  -  CONSTANT  LOAD  OF  FINITE  DURATION 

Statement  of  the  Problem 

In  this  section  the  effect  of  small  damping  on  critical  conditions 
and  estimates  for  dynamic  stability  of  the  models  are  Investigated.  The 
problem  of  constant  load  for  finite  duration  is  considered  here  as  stated 
in  Section  V,  but  small  damping  is  introduced  by  a  dashpot  connected  in 
parallel  with  a  spring. 

Since  damping  is  small,  the  trajectory  of  the  damped  system  could 
be  thought  of  as  a  perturbation  around  the  trajectory  of  the  undamped  one. 
This  is  the  key  assumption  which  helps  find  the  estimates  for  the  damped 
system  as  a  sum  of  the  corresponding  estimates  for  the  undamped  system 
plus  a  product  of  the  small  damping  coefficient,  p,  times  a  quantity 
depending  only  on  the  model  geometry. 

General  Approach 

Damping  forces  F  are  introduced  by  a  dashpot  and  are  equal  to 

F  -  HL 

where  p  is  the  damping  coefficient  and  L  the  trajectory  of  the  dashpot. 

A  nondimensionalization,  p,  of  the  damping  coefficient  p  could  be 
obtained  by 

1 

*  =  XTT  n- 

where  k  is  the  elastic  constant  of  the  spring  and  X  a  normalizing  coefficient 
having  dimension  of  time. 


92 


where  C  is  a  constant  Eq.  (7.2)  characterizes  the  energy  balance  during 
the  action  of  the  load,  p,  and  Eq.  (7.3)  after  the  release  of  P. 

Moreover,  the  continuity  and  both  kinetic  and  dissipated  energies 
at  the  release  of  the  load  yields  the  constant  coefficient  C,  equal  to 


With  this  value  of  C,  Eq.  (7.3),  the  balance  of  energy  after  the  action 
of  the  load  P  becomes 

T°  +  uj  +  5  +  [u*  -  U°]  I  '0  for  t o  <  t  (7.4) 

i  To 

Since  a  system  is  called  dynamically  stable  if  it  is  not  allowed 
to  reach  the  zero- load  unstable  equilibrium  point,  L°  ,  with  the  least 
potential  energy,  the  total  energy  of  the  system,  E°«  0°  +  T°,  with  zero 
velocity  (Kinetic  Energy)  at  L^,  must  be  less  than  the  potential  energy  of 
the  system  at  that  point. 
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Consequently  the  balance  of  energy  for  t  >  tq  expressed  by  Eq. 
(7.4)  gives  the  critical  condition 


(7.5) 


Furthermore  the  nondimcnsionaliecd  damping  coefficient  |I  is  very 
small  (p  «  1)  and  the  trajectory  of  the  damped  system  could  be  thought 
of  as  a  perturbation  around  the  undamped  one.  The  critical  displacement, 

p 

L£r>  the  displacement  which  guarantees  stability  as  long  as  the  applied 
force  is  released  before  the  system  reaches  it,  may  be  expanded  in  a 
Taylor's  series  as. 


P  P  -  P  ? 

L  “  _L r  +  p  ,L  +o(M),t 
cr  0  cr  K  1  cr  ^ 


(7.6) 


where  is  the  critical  displacement  for  the  same  problem  but  without 

the  influence  of  damping. 

Use  of  Eq.  (7.1)  into  Eq.  (7.5)  yields 


U° 


T  L 


u 


-»  J  I, 


cr 


Invoking  the  expression  for  L*  from  Eq.  (7.6),  the  quantity 


K  -  *?] 


may  be  expanded  in  Taylor's  series  also. 


cr 


K  •  "?]  I  •  Or  ‘  e?} 


.1 

- 

+  ^Z“ 


cr 


L  i  -1 
o  cr 


re  -  qi 


('''‘cr)*0^2, 


A 


cr 


(7.7) 


(7.8) 


where  j  stands  for  the  degrees  of  freedom  of  the  model. 


Introducing  Eq.  (7.8)  into  Eq.  (7.7)  and  separating  terns  of  the 
same  order  of  magnitude,  yields  two  equations 


•  *?]  I.. 


(7.9) 


and  for  one  degree  of  freedom  trajectory 


iC  -  - 


L 

r u  o 

0  o*f$P)do* 

«r  y'~M'  - 

!"!  '  “?]>  l!. 


(7.10) 


After  evaluating!.  ,  the  critical  time  t  for  the  release  of  the 
ter  °cr 

load  P  is  found  by  the  balance  of  energy  equation  (Eq.  7.2). 

It  is  simply  mentioned  that  Eq.  (7.9)  gives  the  displacement  estimate 

LP  for  the  undamped  system  and  it  is  the  same  as  that  given  in  Section  V. 
o  cr 

On  each  model  Eqs.  (7.9)  and  (7.10)  will  he  applied  and  the  correction 

of  the  displacement,  ,L  (1*  =  i) .  will  he  found.  Furthermore,  corrections 

1  cr 

for  the  critical  time,  t  *  will  be  computed  also. 

°cr 

Model  A  -  Geometrically  imperfect  Model 

Consider  model  A  exhibited  in  Section  II  hut  under  the  influence  of 
small  damping.  Damping  is  introduced  by  a  dashpot  connected  with  the  spring 
as  shown  on  Fig.  7-1. 


Fig.  7-1  Geometry  of  Model  A  with  Damping 


If  i*  is  the  damping  coefficient,  the  nondimcnsiona Listed  dissipated 
energy,  D,  is  expressed  by 


o 

It  is  assumed  that  the  velocity  of  the  dashpot 
joint  B,  and  that  <p  is  the  trajectory  from  0o  to  6. 
Applying  the  nondiraensionalization 


(7.11) 

is  equal  to  that  of 


T 


£ 


and 


(7.12) 


and  recalling  that  the  unstable  position,  on  the  zero  load  poten¬ 
tial  for  the  undamped  system  is  given  by  9  =  —  ,  the  Eqs.  (7.9)  and  (7.10) 
yield 


9  »  arc  cos  (cos  6 

»  l  c 

cr 


(/ 2  -  Jl  +  sin  eQ)% 

P  J 


(7.13) 


and 


sin 


ep 

o  cr 


(7.14) 


The  evaluation  of 

I  = 


Recalling  the  energy  equation  for  the  undamped  system,  Eqs.  (3.1)  and 

(5.2) 

.2  , _  _ _  2 

09  +  [vl4*inoq)  -  ./I4sin0o  ]  -  p  (cosOQ  "  cos  0)  =  0  Go<ocp<o9  (7.16) 


the  integral  below  is  needed  in  Eq.  (7.14) 

•  2i'Z  <£\ 
in  l  .  - 


n 

r 2  9 
v 

o 


sin 


~2j  d<P 


(7.15) 


and 


oq?2-*-C'/ 1+sinlp  -  yi+sinBo]2  -  p(i O0C/  V  <  -y  (7.17) 


Hence  the  integral  I  becomes 


0  -cos  cr>)  -  [/l+s in  p  -  /l+s  inO  ]“sit\2(y  -°-y\  d  q>  + 
o  o  o  o  \4  2 1  oT 


-T,2  •  2  n  o* 


F  vp (cos  0  -cos  0  PJ-  [/H-sin  cp  -  J l+s  inO  ]  sin  (7  “  *V*  )  d  cp 
J  oocr  o  o  4  2  oT 

«P 


The  critical  time  t  may  be  found  by  invoking  Eq.  (7.2)  for 

cr 

0  <  r  <  T  .  Thus 
o 


j  ,  f  ,  _ . djL--  - _  — 

°cr  0  Vp(cos0  -cos0)  -[./l+sinO-  x/l+sinG"  ]2- p.  f^5rsin2(y-  —Help 
o  o  o  oJ  0  2  Jt> 


(7  18) 


(7.19) 


Expanding  t  in  power  series  of  ji  we  have 
cr 


t  T  +  u,  .  t  .  _  ,-2v 

o=0o  ^  1  o  -*  0  (n.  ) 
cr  cr  '*»* 


(7.20) 


where 


8 

t  ,-° cr  de 

OO  =1  — — - .  -  -  - ==rrrr  - - J.. -  -----V, 

cr  8q  N/p(cos  &o~  cosO)  -  [A  +■  sin0  -  *Ji  +  sinO^] 

is  the  critical  time  for  the  undamped  system,  and 


(7.21) 


^Toi  1  po^cr 


T  _  __o|  ;  i  |>c  -- _ - _ ^  + 

*  °cr  ^  I*  H«”  0  2  6  [p(cos6  -cos6)-(«/l+sin0  -  */l+sin9  ) 2]  ^ 


ire  0^sin2(f"T)do'p}de 


p 

Note  thatjQ^  is  a  function  of  parameters  related  to  the  undamped 

system  only.  The  same  holds  true  for  ^ 

The  governing  equations  for  finding  critical  conditions  in  the 

presence  of  small  damping  ((I  «  1)  arc  Kqs.  (7.13),  (7.14),  (7.21)  and  (7.22). 

Note  that  Eq.  (7.18)  is  employed  in  Eq.  (7.14).  These  four  equations  relate 

the  given  small  damping  coefficient  y,  the  applied  load  p,  and  the  time  parame¬ 
ters  0P  and  6P  .  A  critical  condition  is  expressed  in  terms  of  a  load 
o  cr  1  cr 

level  p  and  the  corresponding  time  Tq  =0t0  •  Thus,  a  critical  con- 

cr  cr  cr 

dition  may  be  found  by  posing  the  problem  as  follows:  for  a  given  small 
damping  coefficient  u  and  load  level  p,  find  (through  the  simultaneous  solu¬ 
tion  of  the  four  governing  equations)  the  corresponding  critical  time  para- 

P  P 

meters,  qT q  and  q  ,  and  position  parameters,  8^  and  0  .  Note  that 

cr  cr  °  c  cr 

the  range  of  p-values  (assigned)  must  be  greater  than  dynamic  critical  load 

for  the  case  of  a  suddenly  applied  constant  load  of  infinite  duration,  with¬ 
out  damping.  The  computational  procedure  involves  the  following  steps:  (a) 
assign  a  p-value  and  compute  J) ^  from  Eq.  (7.13),  (b)  Employ  Eq.  (7.18)  in 

p 

Eq.  (7.14)  and  solve  for  j®cr»  (c)  From  Eq.  (7.21)  solve  for  T  ,  and 

cr 

finally  (d)  Employ  Eq.  (7.22)  and  solve  for  .  A  computer  program  is 

cr 

written  to  accomplish  the  solution  and  numerical  results  are  generated  for  three 
values  of  the  imperfection  parameter  0  (0.005,  0.010,  and  0.020).  These 
results  are  presented  on  Table  7.1. 
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Table  7 


Critical  Conditions  for  Count an!  Load  of  Finite  Duration  in 
the  Presence  of  |iani|iin};  (Model  A) 


p 

o  o 

er 

lTo 

cr 

0.45 

31.25 

29.61 

0.50 

19.55 

11.15 

1.00 

8.97 

1.99 

3.00 

3.42 

.55 

5.00 

2.63 

0.35 

10.00 

2.11 

0.20 

50.00 

0.63 

0.06 

0.41 

33.23 

78.53 

0.45 

21.14 

12.85 

0.50 

12.32 

7.30 

1.00 

7.78 

1.64 

5.00 

2.00 

.31 

10.00 

1.24 

.18 

50.00 

0.42 

.07 

0.44 

13.69 

7.47 

0.55 

9.98 

3.62 

0.65 

6.87 

2.53 

1.00 

5.85 

1.30 

5.00 

1.91 

.27 

10.00 

0.96 

.16 

50.00 

0.39 

.06 

0.005 


0.010 


0.020 


99 


Note  that,  since  a  critical  condition  corresponds  to  a  set  of 

p.  t  values,  a  small  damping  coefficient  p.  lias  a  s l nbi  1  is Ing  effect. 
r  o 

cr 

This  effect,  though,  is  very  small.  For  instance,  at  the  high  values 

of  the  load  p  (say  for  90a  0.00S,  p  3  10)  the  corresponding  value  for 

T  (if  £  *  0.04)  is  2.110+  0.008  =  2.118.  Remember  that  the  p  -  t 
°cr  cr 

curve  for  the  undamped  system  (see  Fig.  r>)  is  very  steep  at  the  high  p- 

value  and  virtually  flat  at  the  low  values  of  p.  On  the  other  hand, 

when  p  *  0.45  (a  value  close  to  p  3  0.432)  the  corresponding  critical 

cr® 

time  is  t  ■  31.25  +  1.184  3  32.434.  Since  the  curve  is  very  flat  at  this 
o 
cr 

load  p- value ,  one  may  conclude  that  the  effect  of  small  damping  Is  virtually 
negligible. 

Model  B 

Consider  model  IS  discussed  in  Section  II,  hut  under  the  influence 
of  small  damping.  Damping  is  introduced  by  a  dashpot  connected  with  the 
spring  on  Fig.  7-2. 


P 


Fig.  7-2.  Geometry  of  Model  B  with  damping. 
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If  is  the  damping  coefficient,  the  nond linens ionalized 
dissipated  energy,  5,  is  expressed  by 


1  ,,  2 
2  to 


^  J'  cpdcp 
o 


(7.23) 


It  is  assumed  that  the  velocity  of  the  dashpot  is  equal  to  that 
at  Joint  C  and  cp  is  the  trajectory  from  0  to  6 . 

Applying  the  nondimens ionalization 


a  /k  -  2a  .  .  d 

a"d  ( >  ■» 7 


(7.24) 


and  recalling  that  the  unstable  position,  Lu>  for  the  "zero  load"  potential 
of  the  undamped  system  is  given  by  0°  " ,  then  Eqs.  (7.9)  and  (7.10)  yield 


2p  (1  -  cos  0  +  e  sin  0  )  =  1 

o  cr  o  cr' 


(7.25) 


4. 

r  . 


0  =  _ 0oT  O  T _ 

1  cr  2p  (sin  9  +  t  cos  0  ) 

o  cr  o  cr 


(7.26) 


Thus  for  the  evaluation  of  the  integral 

n 

2 

i  =  r  <d  d  9 
V  ° 


(7.27) 


the  balance  of  energy  expressions  for  the  undamped  system  Eqs.  (5.1)  and  (5.4) 
are  recalled .  Hence 


•  2  2 

cp  +  sincp-2p(l-  cos  cp  +  e  sinocp)  = 


0  for  (X  cp<  0 

ox  o  cr 


(7.20) 


cp2  +  sin  cp  -  2p  (1  -  cos  0  +  c  sin  0  )  *  0  for  0<  cp<^  (7.29) 

o  ox  o  cr  o  cr  o  croT  2 
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Integral  I  becomes 


f  “  & 
0 


(1  -  cos  9  +  c  sin  9)  -  sin  9  d  9  + 
o  or  or  o 


n 

2 


l  > 


(1  •  cos  8  +  e  sin  6  )  -  sin  ip  d  <P 

'  /«  or  rt  or'  q'  q  T 


o  cr 


o  cr 


o  cr 


The  critical  time,  t  ,  may  be  found  by  invoking  Eq.  (7. 

cr 


cr 


d8 


t  **  r  fr  n  w  2  ~  *  0  q 

°cr  0«/2p  (1  -  cos  0+esin8)-sin0-piJ  9d9 
v  0 

Expanding  r  in  power  series  of  p  (p  «  1)  one  obtains 
cr 

2 

T  *  T  +  p.  -T  +  O  (p  ) 

o  00  ^  1  o 

cr  cr  cr 


where 


ep 

o  cr 


T 

o  o 


cr 


-/ 


d0 


°  72p  (  i 


cos  8  +  c  sin  0)  -  sin  0 


is  the  critical  time  for  the  undamped  system,  and 

0P  ,8 

^0  lj?" _ Jqo*  do* _ 

cr  dpi  p=0  2  0  t2p(l  -  cos  8  +e  sin0)  -  sin20]^2 


lTo 


d0 


.•P 

1  cr 


1 


2p  (1  -  cos  0  +  c  sin  0  )  -  sin  0 

r  '  o  cr  n  rr  o 


o  cr 


o  cr 


Note  that  is  a  function  of  parameters  related  to  the 

1  cr 

system  parameters  only.  The  same  holds  true  for  ^ 


(7.30) 

2) .  Thus 

(7.31) 

(7.32) 

(7.33) 

(7.34) 

undamped 


The  governing  equations  for  finding  critical  conditions  in  the  presence 
of  small  damping  (u  «  1)  arc  Eqs .  (7.25),  (7.26),  (7.33)  and  (7.34). 

Note  that  Eq.  (7.30)  is  employed  in  Eq.  (7.26).  These  four  equations  relate 
the  given  small  damping  coefficient  p,,  the  applied  load  p,  the  time  para¬ 
meters  T  and  T  ,  and  the  position  parameters  0P  and  ,8P  .  A 
o  o  1  o  o  cr  1  cr 

cr  cr 

critical  condition  is  expressed  in  terms  of  a  load  level  p  and  the  corre¬ 
sponding  time  T  =  t  +  m  t  .  Thus,  a  critical  condition  may  be  found 
ooo^lo 
cr  cr  cr 

by  posing  the  problem  as  follows:  for  a  given  small  damping  coefficient  u 
and  load  level  p,  find  (through  the  simultaneous  solution  of  the  four 

governing  equations)  the  corresponding  critical  time  parameters,  T  and 

P  P  cr 

r  ,  and  position  parameters,  G  and  ,0  .  Note  that  the  range  of 

1  o  o  cr  1  cr  ° 

cr 

p-values  (assigned)  must  be  greater  than  dynamic  critical  load  for  the  case 
of  a  suddenly  applied  constant  load  of  infinite  duration,  without  damping. 

The  computational  procedure  involves  the  following  steps:  (a)  assign  a 
p- value  and  compute  ^ ^  from  Eq.  (7.25),  (b)  employ  Eq.  (7.30)  in  Eq.  (7.26) 

p 

and  solve  for  ^0  ,  (c)  from  Eq.  (7.33)  solve  for  t  ,  and  finally 

°  °cr 

(d)  employ  Eq.  (7.34)  and  solve  for  t 

cr 

A  computer  program  is  written  to  accomplish  the  solution  and  numerical 
results  are  generated  for  three  values  of  the  eccentricity  parameter  e 
(0.005,  0.010,  and  0.020).  These  results  are  presented  in  Table  7.2. 

Note  that,  since  a  critical  condition  corresponds  to  a  set  of  p,  t 

o 

cr 

values,  a  small  damping  coefficient  yi  has  a  stabilizing  effect.  This 

effect,  though,  is  very  small.  For  instance,  at  the  high  values  of  the 

load  p  (say  for  c  *  0.005,  p  *  10)  the  corresponding  value  for  T  (if 

cr 

^  *  0.04)  is  0.16  +  0.0064  ■  0.1664.  Remember  that  the  p-T  curve  for 


the  undamped  system  (see  Fig.  5.3)  is  very  sleep  at  the  high  p-value  and 

virtually  flat  at  the  low  values  of  p.  On  the  other  hand,  when  p  ■  0.95 

(  a  value  close  to  p  *  0.948),  the  corresponding  critical  time  is 

cr 

CD 

t  •  1.93  +  0.055  *  1.985.  Since  the  curve  is  very  flat  at  this  load 

o 

cr 

p-value,  one  may  conclude  that  the  effect  of  small  damping  is  virtually 


negligible 


Table  7.2.  Critical  Conditions  for  Constant  Load  of  Finite  Duration  in 
the  Presence  of  Damping  (Model  1$). 


e 

P 

T 

o  o 

cr 

lTo 

cr 

0.005 

0.95 

1.93 

1.39 

1.20 

1.28 

1.04 

1.50 

0.75 

0.80 

5.00 

0.23 

0.27 

10.00 

0.16 

0.16 

100.00 

0.029 

0.04 

0.01 

1.15 

1.22 

10.04 

1.70 

0.65 

2.06 

5.00 

0.27 

0.33 

10.00 

0.18 

0.14 

20.00 

0.11 

0.07 

100.00 

0.024 

0.06 

0.02 

1.10 

0.78 

13.96 

1.50 

0.46 

2.90 

5.00 

0.16 

0.33 

10.00 

0.092 

0.14 

30.00 

0.048 

0.04 

75.00 

0.022 

0.02 
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Model  C 


In  this  article,  the  effect  of  small  damping  on  a  two-degree  of  freedom 
model  (model  C)  is  investigated.  Damping  is  introduced  through  a  dashpot 
connected  with  the  extensional  spring  as  is  indicated  on  Fig.  4.2. 

The  damping  coefficient  is  equal  to  p. 


the  nondimens ionalized  damping  coefficient  and  fj  the 
nond imens ionn lized  dissipated  energy,  then 


D 


;  j  £  d“ 


(7.37) 


and  with  the  right  hand  sides  of  the  Eqs.  (7.3b)  and  (7.36),  the  non¬ 
dimens  tonal  ize  dissipated  energy  becomes 


5  *  £  £  (2rr  +  6ss)  (2rdr  +  6sds)  (7.38) 

where  §  stands  for  the  line  integral  along  the  trajectory  of  the  system 
In  the  r-s  plane.  Furthermore,  the  lowest  potential  unstable  equilibrium 
position  of  the  "zero  load"  potential  is  the  saddle  point  with 

“  f  (A  -  1)  (7.39) 


Recalling  that  for  the  undamped  system  the  trajectory,  which  must  be  followed 
for  maximizing  the  added  energy  when  p  is  acting,  lies  on  the  r  axis,  and 
recalling  Eq.  (7.9),  it  is  found  that 


P 


r 

o  cr 


.2 

4p 


(A  -  1)  +  V/A 


(7.40) 


However,  the  trajectory  that  the  system  follows  from  the  release  of  the 
load  until  it  reaches  the  unstable  saddle  point  (r=>/A--|,s»  +  (\/A%^)^) 
is  not  known.  Since  Eq.  (7.38)  gives  the  dissipated  energy  during  this 
period  of  time  as 


D  ■  p 


7a-§ 
t  p 

r 

o  cr 


(2r  +  6ss  )^  r  dr 


a  lower  bound  for  the  absolute  value  of  D  may  be  found,  if  one  considers 


the  path,  e  =  0  and  j  '  r  ft.  Note  that  r  ■  ft  Is  the  starting 

point  and  r  ■  ft  -  —  is  the  r-coordinate  of  the  saddle  point  on  the  "  zero 
load"  potential. 

Following  the  same  steps  as  in  the  brachistochrone  problem  (see  Append  Lx 
A),  it  is  proven  that  the  symmetric  path  s  0  is  the  minimizing  path  for 


5  min.  Then  from  Fq.  (7.10) 


1  cr  2p  •> 

ft 


(7.41) 


where  r£r  is  the  parameter  shown  in  the  equation 


rP  =  rP  +  (I  .rp  +  o  (|I2) 
cr  o  cr  1  cr 


(7.42) 


Furthermore,  recalling  that  the  total  energy  for  the  undamped  system  along 
the  r  axis  (  s  «  0)  is  expressed  by 


\  (r)2  +  (r  -  ft )2  +  |  (A  «  r2)2  -  2P  (ft  -  r)  *  0  , 


(7.43) 


?  =  -  J&p  (ft  -  r)  -  (A  -  r2)2  -  2  (r  -  ft)‘ 


(7.44) 


Hence,  Eq.  (7.41)  becomes 
P 


1rP^  “  ’  p  f  f2  v4p  ^  ■  r>  "  “  r2y  *  2  (r  -  ft) 


V  dr 


■ft  3 

2%/?'2  2 


-iS  P  r  >  <^A  -  orcr>  •  (A  -  r  >  "  K'-ft)' 


(7.45) 


In  addition,  the  critical  time  tq  may  be  found  by  invoking  Eq.  (7.2),  which 
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expresses  the  balance  of  energy  for  0  <  t 


C^A-r)  -  (A-r2)2  -  2  (r -^A)  -  *  J^x2  dx 

Jr 


and  expanding  tq  in  power  series  of  6t  «  1),  one  may  write 


T  -  T  +  p  T  +  0(0 
O  O  O  1  o 

cr  cr  cr 


(7.47] 


where 


t  -  r  |  =  ’  r  - 

j  n  o  1  J  _ 


°  °cr  °cr  p-0  0  /4p  <JK-r)  -  (A-r2)2  -  2  (r  -,/ft)2 


(7.48] 


°cr  ,  °rcr  JR _ 

**  p*0  i  [4p  (JR  -  r)  -  (A  -  r2)2  -  2  (r  - -A)]3/2 


rP  8J  x2  /4P  C^-*>  “  <A-x2)2  -  2  (x-^/Jv)2  d: 


/4p  6/A  -  r  )  -  (A  -  r  )Z  -  2  (  r  ->/R) 
'  o  cr  '  o  cr'  'o  cr  '  ' 


(7.49) 


Note  that  .r  is  a  function  of  the  undamped  system  only.  The  same  holds 
i  cr 

true  for  .t 

°cr 

The  governing  equations  for  finding  critical  conditions  in  the  presence 

of  small  damping  (p<<  1)  are  Eqs.  (7.40),  (7.45),  (7.48)  and  (7.49). 

These  four  equations  relate  the  given  small  damping  coefficient  m  the 

applied  load  p,  the  time  parameters  T  and  jT  ,  and  the  position 

p  p  °cr  cr 

parameters  r  _  and  ,r  .  a  critical  condition  is  expressed  in  terms  of  a 
o  cr  l  cr 

load  level  p  and  the  corresponding  time  t  =  t  +  ^.t  .  Thus,  a 

°cr  °  °cr  1  cr 


A 


9 


critical  condition  may  be  found  by  posing  the  problem  as  follows:  for  a 

given  small  damping  coefficient  ^  and  load  level  p,  find  (through  the 

simultaneous  solution  of  the  four  governing  equations)  the  corresponding 

critical  time  parameters,  t  and  t  ,  and  position  parameters,  rP 
,,  o  °cr  1  °cr  o  cr 

and  ^r„r*  Note  that  the  range  of  p-valucs  (assigned)  must  be  greater  than 
dynamic  critical  load  for  the  case  of  a  suddenly  applied  constant  load 


of  infinite  duration,  without  damping.  The  computational  procedure  Involves 

p 

the  following  steps:  (a)  assign  a  p-valuo  and  compute  rcr  from  Eq.  (7.40), 

P 

(b)  employ  Eq.  (7.45)  and  solve  for  ^r^^  (c)  From  Eq.  (7.48)  solve  for 

t  ,  and  finally  (d)  employ  Eq.  (7.49)  and  solve  for  T 
o  o__  1  o 

er  cr 

A  computer  program  is  written  to  accomplish  the  solution  and  numerical 


results  are  generated  for  three  values  of  the  parameter  A  (A  =  5.0,  6.0,  8.0). 
These  results  are  presented  on  Table  7.3.  Note  that,  since  a  critical 


condition  corresponds  to  a  set  of  p,  T  values,  a  small  damping  coefficient 
-  cr 

p.  has  a  stabilizing  effect.  This  effect,  though,  is  very  small,  For 


instance,  at  the  high  values  of  the  load  p  (say  for  A=5,  p=10)  the  corresponding 

value  for  Tq  (if  a -0.04)  is  0.32  4-  0.029  -  0.349.  Remember  that  the 
cr 

P-tq  curve  for  the  undamped  system  (see  Fig.  5.5)  is  very  steep  at  the 
cr 

high  p-value  and  virtually  flat  at  the  low  values  of  p.  On  the  other  hand, 

when  p=4.15(a  value  close  to  pcf  “3.71)  the  corresponding  critical  time 

is  r  =  1.32  +  3.18  *  4.50.  Since  the  curve  is  very  flat  at  this  load 
o 

cr 

p-value,  one  may  conclude  that  the  effect  of  small  damping  is  virtually 
negligible. 
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Tabic  7,3.  Critical  Conditions  for  Constant  l.oad  of  Finite  Duration 
in  the  Presence  of  Dampiii}\  (Model  C) 


p 

T 

°  °cr 

lTo 

cr 

4.15 

1.32 

79.66 

4.20 

1.29 

33.11 

4.30 

1.13 

16.50 

5.00 

0.78 

4.12 

8.00 

0.39 

1.10 

10.00 

0.32 

0.74 

30.00 

0.11 

0.10 

70.00 

0.41 

0.03 

100.00 

0.29 

0.02 

5.25 

1.70 

131.11 

5.50 

1.52 

5.12 

7.00 

0.60 

2.94 

10.00 

0.39 

1.57 

12.00 

0.34 

0.78 

20.00 

0.18 

0.31 

30.00 

0.13 

0.16 

50.00 

0.078 

0.06 

80.00 

0.042 

0.03 

7.75 

0.71 

56.31 

8.50 

0.65 

6.12 

9.00 

0.55 

4.34 

15.00 

0.23 

0.99 

20.00 

0.21 

0.68 

30.00 

0.14 

0.34 

50.00 

0.089 

0.11 

80.00 

0.052 

0.05 

150.00 

0.018 

0.01 
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APPENDIX  A.  THE  BRACHISTOCHRONE  PROBLEM 

Constant  Load  of  Finite  Duration-Model  C. 

The  problem  of  dynamic  stability  of  model  C  (Snap-through  Model) 
subjected  to  constant  load  of  Finite  Duration  and  discussed  in  Section  V, 
reduces  to  finding  the  brachistochrone  path  starting  from  point  (r,  s)  = 

(t/K,  0)  and  reaching  the  vertical  r  =  rcr  *JR  -  (A-l)  (A.  1) 

where  p  is  the  nondlmensionalized  load  parameter  and  A  is  a  geometric 
parameter  of  the  model.  Since,  from  physical  considerations,  the 

motion,  in  the  interval  a/A-3  <  r  <  Jh  follows  the  steepest  descent  path, 
which  is  the  symmetric  path  (s=0) ,  the  brachistochrone  problem  will  be 
restated  as  follows:  find  the  path  that  requires  the  smallest  time,  for 
the  system  to  reach  positions  characterized  by  r  *  r  and  starting  from 
(r  =■  JK- 3,  s*0)  with  s#=0. 

The  total  Potential  Energy  of  the  system  is  given  by  Eq.  (2.13)  as 

u£  *  (r2+9s2-2^r+A)  +  \  (A-r2-3s2)2-  (A.  2) 

Recalling  Eq.  (5.1). 

UP  +  TP*  0  for  0  <  t  s  Tq  (A.  3) 

where  TP=  (1+s'2)  (^)2  (A. 4) 

and  <>’-&• 

then  the  brachistochrone  problem  reduces  to  extremizing  the  following  integral 

rcr  .  r~T  rcr 

I-J  (*  “  J  *  ('■*■•'  (A.5) 

Jfc 3  '2ut 

by  considering  paths  which  start  it  r  >  */A-3  and  s  ■  0  with  s'-  0  and  reaching 
positions  characterized  by  r  -  rcj.  (constant). 
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The  Euler-Lagrange  equation  and  the  associated  boundary  conditions  for 


this  brachistochrone  problem,  Eq.  (A. 5),  are  (see  reference  [26] 


9X  _  _d_  =  n 

?)s  dr 


(A. 6) 


dX 

^  s 1 


0  -  s' 


0 


Note  also  that  s  =  s'  =  0  r  *  . 

Since  X  =  -  V — — r-  ,  Eq.  (A. 6)  yields 
-2lT 


s"  +  (s'  +  a 1 3) r+9ss '  VA-(A-r2-3s2)  (r+3ss Hp 

'  '  — p 


-u; 


.  (u,,yw!MVi.0 


(A.  7) 


u: 


One  solution  to  Eq.  (A. 7),  that  satisfies  the  boundary  conditions,  is  given 

by  the  path  s=0.  However,  there  is  no  guarantee  that  path  s^O  is  the  unique 

extremal  path  of  the  variational  problem.  Suppose  that  the  solution  to  Eq. 

(A. 6)  bifurcates  at  a  point  (r  *  r.  s  -  0)  where  r  <  r,  <  Jh-3.  Obviously 

d  cr  —  D  — 

s'«  0  at  that  point. 

However, considering  a  closed  domain  D  =  [r,s,s'],  such  that  s  <  6  and 
s'  <  fi  where  6  is  a  very  small  positive  number  (6«1)  with  the  point  (r^  , 

— p 

0,  0)  lying  in  the  interior  of  D  and  U^,  <-  w  (w  >  0  any  proper  positive  number) 
everywhere  on  D, the  function 

f(r,s,s')  «- (s *+s ' 3)-qss L?/E-L^2-l*2)  C^3ss  1  Hp 


-(1+8 


•  2.  2  9s-3(A-r2-3s2) 

'  — D 


u: 


(A. 8) 
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is  continuous  and  satisfies  the  following  Lipschltz  condition  on  D. 


|f(r,s,s')  -  f(r,s,s,)l<L1|s,-s'j+L2js-sl 

Since  higher  order  terms  in  s  and  s'  are  cancelled  out  and  and  L2  are 
equal  to 

L.  »  max  )r+P)  =  maxQ^r.s) 

UT 

and  L_  =  max  |  — — =  max  Q?(r,s) 

The  existence  of  and  L2  18  guaranteed  by  the  continuity  of  the  functions 
Qj  and  fl2  on  D.  See  reference  [273. 

II 

Recalling  the  Uniqueness  theorem  (Picard-Lindelof)  for  a  differential 

equation  of  second  order,  see  reference  [28],  the  initial  value  problan, 

composed  by  the  differential  Eq.  (A. 7)  and  the  boundary  conditions  s|  = 

rb 

s' |  =  0,  used  as  initial  conditions  in  the  mathematical  sense,  admits  a 

rb 

unique  continuous  solution  in  the  interior  of  D.  Therefore,  the  solution 

to  Eq.  (A. 7)  which  is  initially  symmetric  does  not  bifurcate  at  any  point  on 
the  r-axis.  Consequently  the  synmetrlc  solution  is  unique. 

P 

Since  the  total  potential  is  symmetric  with  respect  to  the  r-axis,  the 
extremal  path,  s=0,  cannot  be  a  minimax.  Moreover,  around  any  path  there 
always  exists  another  path  yielding  higher  value  for  the  time  of  motion 
[integral  I  expressed  by  Eq.  (A. 5)].  Thus,  the  extremal  path  s=0  cannot  be  a 
maximal  path  either.  Therefore,  s=0  is  the  brachistochrone  path  of  the  problem. 

Note  that  path  s=0  is  not  a  ficticious  path  but  a  real  one,  since  it 
satisfies  the  equations  of  motion. 
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